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ABSTRACT 


ANALYSIS AND TESTING OF PLATES WITH PIEZOELECTRIC 
SENSORS AND ACTUATORS 

Jeffrey S. Bevan 
Old Dominion University, 1997 
Director: Dr. Chuh Mei 

Piezoelectric material inherently possesses coupling between electrostatics and 
structural dynamics. Utilizing linear piezoelectric theory results in an intrinsically 
coupled pair of piezoelectric constitutive equations. One set describes the direct 
piezoelectric effect, where strains produce an electric field, and the other set describes the 
converse effect, where an applied electrical field produces strain. The purpose of this 
study is to compare the finite element analysis and experiments of a thin plate with 
bonded piezoelectric material. 

Since an isotropic plate in combination with a thin piezoelectric layer constitutes a 
laminated composite, the classical laminate plate theory is used in the formulation to 
accommodate generic laminated composite panels with multiple bonded and embedded 
piezoelectric layers. Additionally, the von Karman large deflection plate theory is 
incorporated in the stress-strain relations of the laminate. The formulation results in 
laminate constitutive equations that are amenable to the inclusion of the piezoelectric 
constitutive equations, yielding a fully coupled electrical-structural composite laminate. 

Using the finite element formulation, the governing differential equations of motion 
of a composite laminate with embedded piezoelectric layers are determined. The finite 
element model (FEM) not only considers structural nodal degrees of freedom (d.o.f.) but 
an additional electrical d.o.f. for each piezoelectric layer. 
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Comparison is performed by treating the piezoelectric first as a sensor, and then again 
as an actuator. To assess the piezoelectric layer as a sensor, uniformly distributed 
pressure loads were applied and the corresponding generated voltages were determined 
using both linear and nonlinear finite element analyses. Experiments were carried out by 
applying the same uniform distributed loads and measuring the resulting generated 
voltages and corresponding maximum plate deflections. It is found that a highly 
nonlinear relation exists between maximum deflection and voltage versus pressure 
loading. 

The dynamic sensor was evaluated by comparing the predicted sensor voltage with 
the experimental voltage due to a sinusoidal point excitation. In order to assess 
piezoelectric actuation, a sinusoidal excitation voltage was applied and the center plate 
deflection was measured experimentally and compared to the predicted displacement. 
The plate deflection, as a function of time, was determined using the linear finite element 
analysis. 
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FOREWORD 


The research results contained herein partially fulfill the requirements of NASA 
research grant NAG1-1684, entitled "Experimental and Numerical Analysis of Structural 
Acoustic Control for Interior Noise Reduction." This Masters of Science thesis prepared 
by Jeffrey S. Bevan under the guidance of Professor Chuh Mei of Old Dominion 
University, Aerospace Engineering Department, constitutes the research results contained 
herein. The report presents a coupled electrical-structural finite element formulation, 
finite element analysis, and experimental results of a panel with a bonded piezoelectric 
sensor and actuator. The Aerospace Engineering Department, Old Dominion University 
and Langley Research Center Structural Acoustic Branch both provided computational 
and experimental facilities required to complete the research study. Mr. Travis L. Turner 
of Langley Research Center Structural Acoustic Branch was the technical monitor. 
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CHAPTER I 


INTRODUCTION 
1.1 Historical Perspective 

The piezoelectric effect was first identified in 1880 by Pierre and Jaques Curie [1] and 
has remained an active research topic ever since. The Curie's discovered the direct 
piezoelectric effect by placing a weight upon a crystal and observing that a charge 
proportional to the weight was generated. Shortly thereafter they confirmed the converse 
piezoelectric effect by observing an induced strain resulted when a voltage was applied to 
the crystal. Hence the term piezoelectricity, meaning pressure electricity, was coined to 
describe this phenomenon. Piezoelectricity remained somewhat of a scientific curiosity 
since the complexity of the coupled electrical and mechanical properties were unknown. 
Thus one of the objectives of the earliest research efforts was to better understand the 
coupled electrical-structural properties of piezoelectric material and to develop accurate 
analytical models to support and direct engineering design applications. 

A portion of the mysterious veil was lifted from piezoelectricity during World War 1 
when Professor Langevin, under the auspices of the French government, set out to 
determine a method to detect submarines [1]. Professor Langevin used piezoelectric 
crystals in a device that, when submerged underwater, generated a voltage when a 
disturbing wave front would impinge upon it. Conversely when the device was 
electrically excited it would vibrate and emit a longitudinal underwater wave. Professor 
Langevin was unable to conclude his research until after the war, however his device was 
the predecessor of today's modem sonar transducer. 

Another early application of piezoelectricity was discrete crystal circuit devices such 
as oscillators and filters. The crystal oscillators were extremely stable and were used 
extensively in military communication equipment. At one point there were in excess of 
30 million crystals in military equipment in one year [1]. The crystal controlled 
oscillators resulted from the research efforts of Cady at Weslyan University [1]. Not only 
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do piezoelectric crystals possess a characteristic stable resonance, they also are extremely 
selective, which is indicated by their high mechanical Q values. This sharpness provided 
the ability to design extremely discriminate filters, resulting in precision circuitry capable 
of separating simultaneous multiplexed conversations over a single pair of wires. It is not 
difficult to realize the important role that piezoelectric crystals had in the development of 
today’s modem telecommunications industry. 

Another milestone in understanding piezoelectric phenomena was contributed by 
Professor Mindlin. Mindlin began ground breaking analysis on waves and vibrations in 
isotropic elastic plates concurrently with high frequency vibration of crystal plates [2]. 
Subsequent work on isotropic bars and plates lead to unprecedented design and 
development of electromechanical filters and discrete time delaying devices. Mindlin's 
pioneering papers on crystal plates may be considered the most significant in modem 
piezoelectric research, since it clarified the complicated coupled piezoelectric 
phenomena, leading the way to improved piezoelectric designs for quartz crystal filters. 
Mindlin's research lead to a sole-supplier contract in 1955 from the U.S. Army Signal 
Corps, a long time sponsor of the research on crystal plate vibrations, resulting in a 
monograph entitled "An Introduction to the Mathematical Theory of Elastic Plates". 
Since the application of quartz filters and other circuit devices such as surface acoustic 
wave devices, piezoelectric materials have found uses in numerous applications such as 
dot-matrix printers, computer keyboards, high-frequency stereo speakers, igniters, 
microphones, accelerometers, and various transducers (force, strain, and pressure), 
however a new and active research area commonly referred to as smart structures has 
become a very popular research topic. 

The concept of smart structures is a relatively new and diverse field. The 
fundamental core of smart structures integrates sensors and actuators to structural 
elements to obtain a state of desirable static and, or dynamic control [3]. The 
development of smart structures has resulted from three recent significant trends [4]. The 
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first is the increased utilization of traditional laminated composite structures. The 
composite material theory incorporates smaller constitutive elements thus permitting the 
inclusion of piezoelectric constitutive relations with relative ease. It is not unreasonable 
to visualize a structural member consisting of multiple sensors, actuators, and signal 
processors. 

The second trend is the consideration of coupling the structural and electrical 
properties by utilizing the off-diagonal terms of the constitutive relations. This practice 
has been utilized by the large deflection laminated composite theory by including the 
coupled bending-extensional laminate stiffness. The third trend is the advanced rate of 
growth within the computer science and electrical engineering disciplines. Hardware 
improvements such as miniaturization and increased computational power permit faster 
and more complex algorithms resulting in greater adaptability of the smart structure. 

Even though Mindlin's work transformed piezoelectricity from a state of scientific 
curiosity to an applied science, his higher-order theory of wave propagation is not directly 
applicable to low frequency structural vibrations of laminated structures. Considering 
that the physical geometry dictates that the structural resonances occur several orders of 
magnitude below the piezoelectric crystal's natural frequency, linear piezoelectric theory 
is exploited where the only electrical and structural interaction arises from the linear 
piezoelectric constitutive relations [5]. 

The governing equations for distributed piezoelectric using linear piezoelectric theory 
combined with the classical laminate plate theory was presented by Lee [6]. Subsequent 
research has exploited linear piezoelectric theory for numerous smart structure 
applications such as self-sensing piezoelectric actuators [7] and modal analysis using 
piezoelectric sensors [8]. The self-sensing piezoelectric actuator results in a truly 
collocated sensoriactuator capable of simultaneously measuring a structure's dynamic 
response while providing a controlling input. The collocation characteristic provides a 
cost benefit by reducing the number of transducers [9]. Furthermore a classical coupled 
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electrical-structural analysis approach was demonstrated by Lai [10] to control panel 
flutter. 

The classical analytical solution method of piezoelectric structural analysis has been 
extremely useful, however their solutions are restricted to relatively simple geometries 
and boundary conditions. Since the finite element method has proved to be a powerful 
and popular technique for the analysis of complicated structures and multi-field problems 
it may be applied to the coupled electrical-structural piezoelectric system. 

Piezoelectric solid finite elements were formulated and applied to transducers and 
oscillators by Allik and Hughes [11]. However since smart structures typically consist of 
thin piezoelectric layers attached to structures that are several orders of magnitude greater 
in thickness, the solid finite element formulation leads to an inherently inefficient process 
in which to model the complete physical structure. Thus Tzou and Tseng formulated a 
new thin piezoelectric solid finite element coupled to shell and plate finite elements [12]. 
The intrinsic parasitic shear effects associated with such finite elements were eliminated 
by the introduction of internal d.o.f.'s within the piezoelectric plate element. 

One such application of the finite element method incorporating piezoelectric sensors 
and actuator was presented by Zhou in analyzing and suppressing panel flutter [13]. The 
flexibility of the finite element method was demonstrated by Zhou by formulating a 
structural, electrical, and thermal coupled analysis resulting in the control and 
suppression of panel flutter of composite laminated panels at elevated temperatures. 
Zhou provided simulation results for aerodynamically induced large deflections with 
multiple embedded piezoceramic actuators for various ply orientations. In addition Zhou 
introduced a novel time domain nonlinear solution method by developing a set of forced 
Duffing equations in reduced modal coordinates. By introducing the modal 
transformation, the number of equations to be solved is greatly reduced thus affording 
great computational savings. 
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1.2 Objective and Outline 

Much research has been conducted using linear piezoelectric theory for control in a 
variety of structural vibration applications. Both classical analysis and the finite element 
method has been utilized. The objective of this research is to compare the fully coupled 
electrical-structural finite element analysis of an isotropic panel with a surface mounted 
piezoelectric patch to experimental test results. 

This thesis is organized as follows. In Chapter I, a historical background and the 
objective of this study is presented. Chapter II introduces the piezoelectric phenomena. 
In addition, an electrical equivalent model and the coupled linear piezoelectric 
constitutive relations are introduced. Chapter III presents the finite element formulation, 
including the fully coupled electrical-structural constitutive relations. An effort has been 
made to develop the most general formulation applicable to laminated composite plates. 
To this end, the von Karman large deflection theory was incorporated, based on the 
results obtained from the experiments conducted. The finite element formulation is then 
modified to accommodate the inclusion of the electrical degree of freedom to satisfy the 
electrical-structural coupled linear piezoelectric constitutive relations. Next, the finite 
element matrices are derived using Hamilton's principle and the subsequent equations of 
motion assembled. Lastly, solution procedures for the static and dynamic sensor and 
dynamic actuator are presented. Chapter IV presents the experimental aspects and the 
test plate clamping fixture design. The fixture design specifications and boundary 
conditions are discussed along with the applied loading procedures. In addition a 
discussion of piezoelectric wafer preparation and mounting is presented. Chapter V 
presents and compares the predicted and experimental results. Comparison of the results 
obtained from test and analysis of a quasi-static sensor subjected to uniformly distributed 
loading is presented, along with results for a dynamic sensor due to mechanical point 
loading and dynamic piezoelectric actuation. Chapter VI provides a discussion of the 
results and conclusions. 
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CHAPTER II 


PIEZOELECTRICITY 
2.1 Introduction 

When a piezoelectric material is subjected to mechanical strains or stresses it gives 
rise to an electric polarization, or it simply generates an electric charge. This 
phenomenon is referred to as the direct piezoelectric effect. Conversely, a piezoelectric 
material will undergo strain when it is electrically polarized, or subjected to an electric 
field. This action is referred to as the converse or reciprocal piezoelectric effect. It is 
important to note that piezoelectric materials are polarized such that elastic deformation 
depends on the sign and magnitude of the applied electric field. That is to say, 
piezoelectric material may undergo either elongation or contraction, simply by reversing 
the polarity of the applied electric field. It is this polarization that differentiates 
piezoelectricity from electrostriction. Electrostriction is a function of the square of the 
electric field, thus sign independent [14]. Since piezoelectric materials exhibit both direct 
and reciprocal effects the same specimen may be implemented as an actuator or sensor, or 
both simultaneously [7,15]. The reciprocal effect facilitates actuation and the direct 
effect favors sensing of structural vibrations. The piezoelectric effect maintains a linear 
relationship between the electrical and mechanical quantities [5]. Thus linear 
piezoelectric theory couples quasi-electrostatic field equations with a dynamic structural 
system. This quasi-electrostatic approximation is valid since the phase velocities of the 
structural vibrations are several orders of magnitude less than the phase velocities of the 
electromagnetic waves [5]. The direct and reciprocal piezoelectric property constitutes an 
inherent electromechanical coupling that is included in the constitutive relations of the 
structural analysis problem considered herein. 

As stated in Section 1.1 piezoelectricity was first discovered in the 1 800's. The initial 
piezoelectric observations utilized naturally occurring crystals such as quartz, tourmaline, 
and Rochelle salt. Piezoceramics, a manufactured polarized ferroelectric material, are 
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commonly used in commercial transducer applications and will be analyzed and tested in 
this study. 

Even though piezoelectricity is a linear effect, a complex relationship exists between 
elastic, mechanical, thermal, and electrical properties. The relationships are shown 
schematically in Figure 2.1 [14]. The comers of each of the triangles are functions 
representing the E electric field, D electric displacement, <J stress, s strain, T 
temperature, and S entropy. The piezoelectric effect is shown by the left-hand side of 
Figure 2.1 indicated by the independent variables d, g, e, and h, which are the 
piezoelectric coefficients, the dielectric permittivity 6 and impermeability fi, the stiffness 
c and compliance s. Linear piezoelectric theory assumes constant entropy therefore it is 
an adiabatic process, thus mechanical strains do not contribute to the thermodynamic 
state of the piezoelectric. 

2.2 Electrically Equivalent Piezoelectric Model 

The piezoceramic used for this work is in the shape of a thin plate with electroplated 
electrodes on the top and bottom surfaces as depicted in Figure 2.2. Physically the 
piezoceramic consists of two conductors separated by a dielectric material. Electric 
circuit theory calls such a device a capacitor. In circuit theory a capacitor is a passive 
device, however the piezoceramic is a polarized ferroelectric material that generates a 
charge proportional to strain as dictated by the direct piezoelectric effect. Thus an 
electrically equivalent model of the piezoelectric material consists of two charge 
generators, a capacitor and a resistor as shown in Figure 2.3 [7]. The charge generators 
q a , qp shown in Figure 2.3 represent the applied charge and self-induced charge, 
respectively. The resistance Rp is the intrinsic electrical resistance of the dielectric which 
in most cases is very large and represents an electrical open circuit condition which 
impedes the flow of electric current. 

In reference to piezoceramic materials, the applied charge results from an externally 
applied voltage, and since the piezoceramic is capacitive, the charge accumulates on the 
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electrodes. By definition electric current is charge in motion, hence the total current 
flowing out of an enclosed volume must equal the loss of charge within the volume. This 
is a statement of conservation of charge and leads to the explanation of why current flows 
in the leads of a capacitor being charged (or discharged) when no current flows between 
the capacitor plates. Since the piezoelectric specimen behaves as a capacitor the current 
flow through a piezoelectric by an applied constant voltage is analogous to the charging 
or discharging of a capacitor [16]. Current flows across the open circuit plates of a 
capacitor since there is an accumulation of charge on the plates. This can be shown by 
examination of statement of conservation of charge 


<jj- da. \p c dV (2.1) 

where J is defined as current density, that is current per meter width, p c is the charge 

density, and da is the differential area with a unit normal vector. 

If the volume enclosing the charge remains constant with respect to time then the time 
derivative may be moved into the volume integral thus 


= ( 2 . 2 ) 

y & 

By applying the divergence theorem, the surface integral is converted into a volume 
integral so that Eq. (2.2) becomes the time-varying equation of continuity and can be 
written as 


dt 

If Gauss's law, V • D = p c is substituted into Eq. (2.3) it becomes 

V- J=— — V-D 
dt 


(2.3) 


(2.4) 
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The electric displacement density D is frequently called flux density, and a more 
complete description may be found in Appendix C. Interchanging differentiation with 
respect to space and time in Eq. (2.4) yields 

v {? +, )-° < i5 > 

Casting Eq. (2.5) in integral form and applying the divergence theorem results in 


<j| — + j)-</a = 0 


( 2 . 6 ) 


Hence Eq. (2.6) indicates that the total current of time-varying fields is f — + , where 


) 


dD . 


dt 


is a displacement current density due to the time rate of change of the electric 


displacement density and J is the current density resulting from the flow of charge. Thus 

when considering a capacitor being charged with a direct current, the time varying current 

density J is zero due to the open circuit condition, but Eq. (2.6) indicates the existence of 

a displacement current density — flowing through the leads of the capacitor being 

dt 


charged, or the piezoceramic. 

2.3 Piezoelectric Constitutive Equations 

Electric enthalpy density H describes the amount of energy stored within the 
piezoelectric material. Given the electrical-structural coupling of the piezoelectric 
material, the electric enthalpy density is the internal coupled strain energy less the stored 
electrostatic energy density [17]. The stored electrostatic energy density is analogous to 
the structural elastic strain energy. A detailed derivation of the electrostatic energy may 
be found in Appendix C. The electric enthalpy density is defined as 

H = U- DE (2.7) 
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where U is the strain energy, D and E are the electric displacement density and electric 
field, respectively whose product represents the electrostatic energy density. The 
enthalpy may be expanded to yield the following relationship 


so that 


H - \ {s}\Qf {*} - {E} T [e){e} -\{EY [e]' {£} 


dH 


G ' j dS; 


( 2 . 8 ) 


(2.9) 


D;=~ 


dH 

dE; 


( 2 . 10 ) 


where cr , s i} are the stress and strain respectively, D j and E j are the electric 
displacement and electric field respectively, [Q] is the stiffness matrix measured at 
constant electric field (short circuit), [e] G is the dielectric permittivity matrix measured at 
constant strain (clamped), [e\ is a matrix of piezoelectric strain constants, and superscript 
T represents matrix transpose. 

Application of Eqs. (2.9) and (2.10) to Eq. (2.8) produces the following coupled 
electrical-structural piezoelectric constitutive equations 


{<7 } = \Qf{e}-[e] T {E) ( 2 - 11 ) 

p} = MM + [em (2-12) 

Due to practical engineering considerations the piezoelectric strain constants \e\ and 
clamped permittivity matrix [ej are not typically available. However, the stress 
constants [d] and the free permittivity matrix [e]^ are readily available and are related to 
[e\ and [e] £ by the following relations [12] 

[e] = [djQf ( 2 . 13 ) 
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(2.14) 


W = W-[dMVl 

Thus the piezoelectric constitutive equations can be expressed as 

M = \Qf{{s}-[d] T {E}) (2.15) 

{L>} = [rf]{ C 7} + [er{£} (2.16) 

Furthermore Eq. (2.15) may be substituted into Eq. (2.16) yielding 

{£>} = [d\Q] E {{s} - [d] T {£}) + [e]" {E} (2.17) 

The following electromagnetic constitutive relation between the electric displacement 
density and the electric field may be used to clarify the physical meaning of Eq. (2.17) 

{D} = [e]{£} (2.18) 

Thus Eq. (2.17) can be written as 

{£}=[/?r[rf]M + {£} (2.19) 

where WY is the free dielectric impermeability defined by \fiY = [e] a , and {E} is an 
externally applied electrical field. Thus the direct piezoelectric effect results in an 
electric field comprised of two components or sources; one self-generated as shown in the 
first term on the right hand side of Eq. (2.19), and the other due to an externally applied 
voltage as shown in the second term of the right hand side of Eq. (2.19). 
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Fig. 2.1 Elastic, thermal, and electrical properties of piezoelectrics 
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Polarization Direction 



Fig. 2.2 Physical description of the piezoelectric element 
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CHAPTER III 


FINITE ELEMENT FORMULATION 
3.1 Introduction 

The derivation of the governing differential equations of motion for a panel with 
embedded piezoelectric layers is introduced in this chapter. The formulation is based on 
the classic laminate plate theory, including the plane stress assumption and the von 
Karman large deflection theory. The variational energy method facilitates the 
formulation of the linear and nonlinear finite element equations of motion in terms of the 
nodal degrees of freedom (d.o.f.) and the fully coupled structural-electrical properties. 

3.2 Displacement Functions 

The panel with piezoelectric layers shown in Figure 3.1 is modeled using the four- 
node modified C* conforming straight-sided rectangular plate element. Each element 
consists of twenty-four structural degrees of freedom. Each node of the element contains 
four bending d.o.f.'s and two membrane d.o.f.'s to represent the transverse, or out of plane 
and membrane displacements, respectively. A piezoelectric element maintains consistent 
structural d.o.f.'s, however an additional electrical d.o.f. is required to satisfy the 
electrical coupling. Thus the nodal displacement vector is augmented by adding an 
electrical (voltage) d.o.f for each piezoelectric layer present within the element. The 
voltage may either be applied to, or generated by the piezoceramic layer or layers. In 
essence, the electrical d.o.f.'s can be treated as structural displacements. 

The element nodal displacement vector consists of the bending and membrane 
displacements and the voltages which can be written as 

M=H w *Y c 3 - 1 ) 

The bending and membrane displacements of Eq. (3.1) represent the nodal displacements 
and are respectively shown as 


15 


M T 

= {w, 

W 2 W 3 

W 4 

™>xl 

™>x2 

W >*4 




™>y2 


W’yA 

W >xy\ 

W’xy2 W ’xy3 ^xy4 

(3.2) 


= {"1 

Li 2 

u 4 

v 2 v 3 

V a) 

(3.3) 


where w represents the transverse deflection, w, x , w, y are the slopes in the x and y 
directions respectively, w, xy is the second order twist derivative, u and v are the x and 

y membrane displacements, and the numerical subscript denotes the node number. The 
electrical d.o.f.'s represented by {w^\ has voltage components corresponding to each 

piezoelectric layer and will be subsequently described in greater detail in Section 3.3. 

The continuous bending or transverse displacement is approximated using a cubic 
polynomial given as 


w = a } + a 2 x + a 3 y + a 4 x 2 + a 5 xy + a 6 y 2 + a 7 x 3 
+a & x 2 y + a 9 xy 2 + a m y 3 + a u x 3 y + a n x 2 y 2 

+a l3 xy 3 + a u x 3 y 2 + a ]5 x 2 y 3 + a ]6 x 3 y 3 (3.4) 

which can be written in compact matrix form as 

w = L# JM ( 3 - 5 ) 

Similarly the continuous in-plane displacements u and v are approximated using bilinear 
polynomials such as 

u~b x + b 2 x + b 3 y + b 4 xy 

v = b 5 +b 6 x + b 7 y + b g xy (3.6) 

which can also be written in compact vector notation as 
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(3.7) 

= L# vJW 

(3.8) 


The generalized coordinates {< 2 } and {b} maintain a spatial relationship to the nodal 
displacements through coordinate transformations and are functions of time only. The 
transformation relationship is given by 

W = (3-9) 

{b} = KJwJ (3.10) 

Appendix A provides a comprehensive derivation of the bending and membrane 
transformation matrices [T b \, \T m ]. Since the electrical degrees-of-freedom vector 

represents the voltage per piezoelectric layer and does not possess any preferred 
geometrical orientation, coordinate transformation is not required. 

The displacement field may be expressed in terms of the nodal d.o.f.'s by substituting 
the respective coordinate transformations into the displacement field approximations, 
thus substituting Eq. (3.9) into Eq. (3.5) and Eq. (3.10) into Eqs. (3.7) and (3.8) yields 

M = L-ff.JMn} = KUH} (3-H) 

{u) = lHj[T m ]{ W J = [BjwJ (3.12) 

M = (3.13) 

where the [£„],[£„], and [5„] matrices are the shape or interpolation function matrices. 

3.2.1 Linear Analysis 

The majority of the work considered herein is based on the assumption of small 
displacements. The strains are therefore comprised of inplane or membrane strains and 
bending curvatures as 

{£•} = {£°}+ z{k} (3.14) 
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which can be expanded as 


( = ) £ y \ + z \ K y 


where } is the membrane strain vector and {a:} is the bending curvature vector. The 
membrane strains and curvatures are functions of the inplane and bending displacement 
respectively and are defined as 


£ y\ = \ 

r%\ { U >y+ V ’i 


(3.16) 


K v \ \ -W. 


K y > = < -W,yy 
K v \ {- 2w ’*y 


(3.17) 


By introducing the approximations for the displacement field as previously defined and 
recalling the generalized coordinate relationship, the strain vector can be expressed as 

IAJ>, 

{ff°}= lH v j, y W = [cj{6} (3.18) 

Similarly the curvature may be expressed in terms of the generalized coordinates as 


W = ‘ ~L H »\’yy > W = [Q]W ( 3 - 19 > 

where the matrices [C m ] and [C b ] result from the differentiation of the shape functions 
with respect to the dependent variables x and y. Thus [C m ] and [C b ] are defined as 
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(3.20) 


[CJ = 


0 1 0 y 0 0 0 0 
0 0 0 0 0 0 1 x 
0 0 1 x 0 1 0 y 


[Q] = H 


0 0 0 2 0 0 6x 2y 
0000020 0 
000020 0 4x 


0 0 6xy 2 y 2 0 6xy^ 2 y 3 6xy > 

2x 6y 0 2x 2 6xy 2x 3 6x 2 y 6x 3 y 

4 y 0 6x 2 8xy 6y 2 1 2x 2 y I2xy 2 1 8x 2 y 2 


(3.21) 


Expressing the generalized coordinates in terms of the nodal d.o.f.'s the strain and 
curvature vectors may be written as 


{e°} = [C m lT m ]{wJ = KR} (3.22) 

W = [C ( Rh} = [^h} (3-23) 

where [B n ] and [B b are the strain interpolation matrices. Thus the strains in Eq. (3.14) 

can be written in terms of nodal d.o.f.'s as 


{4 = fcW + (3.24) 

3.2.2 Large Deflection Analysis 

The von Karman plate theory is a geometrical nonlinear theory that accounts for 
moderately large deflections and small rotations of the mid-surface of the plate. Thus the 
von Karman large deflection strain-displacement relations are defined as 

{s} = {£ q )+z{k} (3.25) 

where 

{*“} ={*;}+&} (3.26) 

Hence {e ° m } is identical to the previously defined membrane strains {«■“}, and \sl } is the 
membrane strains induced by the large transverse deflection. The strain-displacement 
relations for moderately large displacements is defined as 


19 



(3.27) 
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w, x 
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ry >xx 

£ y 

> = i 

v >y 


> + z< 

— W. > 

9 yy 

r X y . 


!*>y +V >x. 


2w ’* w ’y. 


1 
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In future derivations it will become apparent that it is convenient to express the large 
deflection strain in terms of the slope matrix and slope vector as 


{4}=\[e]{8} ( 3 - 28 ) 

where the slope matrix and vector elements are the derivatives of the transverse 
displacement function. Thus the slope matrix and vector are defined respectively as 


[<?]= 


w. 


0 

0 w,. 




w. 


(3.29) 


< 3 . 30 , 

Utilizing the definition of the slope matrix and vector of Eqs. (3.29-30), the strain due to 
large deflections may be expressed in terms of the generalized coordinates. Recognizing 
that the slope vector is the derivative of the bending shape functions, Eq. (3.28) becomes 




W = cm 


(3.31) 


where [C e ] results from the indicated differentiation and is given as 


[C e ] 


0 1 0 2x y 0 3x 2 2 xy 

0 0 1 0 x 2y 0 x 2 


y 2 0 3 x 2 y 2 xy 2 y 2 3 x 2 y 2 Ixy? 3 x 2 y 2 

2 xy 3 y 2 x 3 2 x 2 y Sxy 1 2 x 2 y 3 x 2 y 2 3x 3 y 2 J 


(3.32) 


Employing the coordinate transformation, the membrane strains can be expressed in 
terms of the nodal d.o.f.'s as 
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(3.33) 


k\=\[olc 0 lUWi = \mAM 

Hence the strain-displacement relations described in Eq. (3.24) can be written in terms of 
the nodal d.o.f.’s as 


{*} = (3.34) 

3.3 Electric Field and Electric Displacement Density 

In Section 3.2 the nodal displacement vector included a term described as the electrical 
degrees of freedom. The most general composite piezoelectric element may consist of 
many piezoelectric layers embedded within a laminated composite panel resulting in a 
electrical d.o.f. vector which contains an applied (or measured) voltage corresponding to 
each piezoelectric layer. The electrical d.o.f.'s can therefore be expressed as 

k}={F, v 2 - vj (3.35) 

where np represents the number of piezoelectric layers presents. 

The electrode of the piezoelectric layer establishes an equipotential boundary 
condition and the dielectric permittivity is assumed isotropic. The voltage establishes a 
linear electric displacement density and electric field through the thickness of the 
piezoelectric material. The electric field strength is the negative of the voltage or 
potential gradient, and is defined as 


E = -VF (3.36) 

The piezoelectric material considered for this study is a thin rectangular plate and 
assumed to be isotropic, thus the stress/charge constants simplify to d 3] ~ d 32 . Since the 

electrodes are on the top and bottom of the piezoelectric plate, polarization occurs only in 
the 3-direction. The stress/charge constant d 33 is assumed to be constant throughout the 

thickness which results in an electric field in the 3-direction only. A detailed explanation 
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of the stress/charge relations is provided in Appendix B. The total electric field due to all 
of the piezoelectric layers can therefore be expressed as 


(3-37) 

where the matrix [B^ ] is a diagonal matrix with elements consisting of the reciprocal of 
the thickness of each piezoelectric layer. Thus may be written as 



0 

i 


(3.38) 


where np represents the number of piezoelectric layers present. 

Summarizing the above results, the generalized strain-displacement relations based on 
the small defection assumption can be obtained by combining Eqs. (3.24) and (3.37) as 


{*}= 


z[B b ] [Bj 0 

0 0 -k] 


w » 

w. 


(3.39) 


Similarly, for large deflections, Eq. (3.39) can be modified by including the large 
deflection strain of Eq. (3.34) yielding 

'*1 




E, 


1 


[«J o 

o o -kL 


w J. 


(3.40) 


Since the electric displacement density is assumed to be generated along the polarization 
axis only (3-direction), Eq. (2.17) can be reduced to 

D 3 =ldjQl{s}-E 3 {d}) + e; 3 E 3 (3.41) 

where the appropriate strain may be employed. The stress/charge coefficients {d} are 
expressed as a vector, due to the geometrical assumptions made during the transformation 
of the piezoelectric constants described in Appendix B. 
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3.4 Constitutive Equations 

The finite element equations of motion for the plate element used will be derived 
using classical laminate plate theory (CLPT). The CLPT is used, since in the simplest 
case a piezoceramic patch bonded to an isotropic substrate constitutes a laminate. In the 
most general case, a laminated composite element will consist of a typical lay-up with a 
number of alternating piezoelectric layers. A typical laminated composite with embedded 
piezoelectric layers is shown in Figure 3.2. A typical isotropic panel with symmetrically 
bonded surface piezoceramic patches is shown in Figure 3.3. The CLPT assumes that the 
piezoelectric is perfectly bonded and that each lamina is in a state of plane stress. For the 
thin plate considered, the rotary inertia and transverse shear deformation effects are 
assumed negligible. 

The stress-strain relations of a specially orthotropic composite lamina and a 
piezoceramic layer is [1 8] 


r > 

<*i 


~Qu 

Q\2 

0 ‘ 
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(3.42) 


and 


r > 
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(3.43) 


where the subscripts s and p indicate the structural and piezoelectric lamina respectively. 
The piezoelectric material considered is assumed to be isotropic in the 1- and 2- 
directions therefore dj i~d^2' The polarization axis of the piezoceramic is assumed to be 
such that a positive strain or elongation in the 1- and 2-directions results from a "positive'’ 
applied voltage referenced to the electrode bonded to the plate. 

The stress-strain relations for the k lh layer of a laminated composite is obtained by 
combining Eqs. (3.42) and (3.43) are 
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Sll 012 016 £ x \d x 

*0y > = 012 022 026 < Sy '-ElAdy ' (3-44) 

Jxy. k _0 J6 0 26 0 66 v ./ kj 

where the transformed reduced lamina stiffness matrix [Q\ is developed from the 

transformation of the principle material coordinates with respect to the global 
coordinates, similar transformations exist for the stress, strain, and the stress/charge 
constants. Appendix B provides a comprehensive derivation of the required principal 
material coordinate transformations. 

For a general orthotropic piezoelectric layer, the generated electric displacement 
density along the polarization axis (5-direction) for the k' h layer may be written as 


011 012 016 

~ [^x dy d xy \ k 0i2 022 026 

016 026 066 


-EAd 


^ K: yhk s -"ik 


Eqs. (3.44-45) may be condensed in matrix form as 

^h=[QWl~E 3 AA) (3-46) 

A* = Ml [Q 1 ({*} - A* {d } k ) + ef 3t A* (3 .47) 

where [q \ and {d} k are the lamina stiffness and stress/charge constants respectively for 

the hfh piezoelectric layer and are transformed to the global x,y coordinates. For a 
composite lamina without a piezoelectric layer, set E 3k - {d} k = 0 . 

3.5 Equations of Motion 
3.5.1 Generalized Hamilton’s Principle 

Finite element equations of motion for the laminated composite panel with fully 
coupled electrical-structural properties are derived utilizing the generalized Hamilton's 
principle [19] to obtain 
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(3.48) 


^8{T-U + W e -W m +W)it = 0 

where T and U are the kinetic energy and strain energy of the system, W e is the electrical 
energy, W m is the magnetic energy, and W is the work done due to external forces and 
applied electric field. The magnetic energy is negligible for piezoceramic materials if no 
external magnetic fields are located near the specimen. The kinetic energy of plate 
element is defined as 


T = J“p({w} r {w} + {u} T {u} + {vf {v})d¥- (3.49) 

v- ^ 

where w, u, and v are the transverse and membrane velocity components and p is the 
mass per unit volume, and — is the volume of the element. The potential and electrical 
energies are defined as 


u= l~{s} r {a}dF 

(3.50) 

W e = f-{Ef{D}d¥ 

V- ^ 

(3.51) 


and the work done on the element by external sources is defined as 

w = fM r {F„ }dV + |M r {F s )dS + {w} T {F c } - \Vp ct dS (3.52) 
v- s, s 2 

where {F b } is the body force vector, {/%} the surface traction vector, {F c } is the 
concentrated loading vector, S { is the surface area of the applied traction, S 2 is the surface 
area of the piezoelectric material, V is the voltage applied to the piezoelectric, and p cs is 
the surface charge density generated by the piezoelectric effect. The electrostatic energy 
results from the charging process of the equivalent piezoelectric capacitance, as described 
in Appendix C. In Hamilton's principle, all variations must vanish at the time t~t x and 
t = t 2 . The Hamilton’s variational statement may be written in the most general form as 
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(3.53) 


J[p({c>w} 7 {w} + {Su} T {u} + {<$v} r {v}) 

V- 

- {Ssf {a} + {SEf {D} + {Swf{F b }}lV 


+ J{w} 7 {F ! }d5- \SVp cs dS + {Swf{F c } = 0 

S l S2 

Evaluation of Eq. (3.53) leads to the development of the finite element matrices and the 
elemental equations of motions. 

3.5.2 Resultant Forces and Moments 

The stresses of each individual lamina are not necessarily equal, therefore Eq. (3.44) 
is not directly applicable since the curvatures are typically unknown and are very difficult 
if not impossible to measure experimentally. However the inplane strains and curvatures 
of Eq. (3.44) can be related to the applied forces and moments through the static 
equilibrium conditions thus making Eq. (3.44) more useful [18]. When working with 
laminated composite plates, it is however very convenient to consider the forces and 
moments per unit length. Such forces and moments are commonly referred to as the 
stress resultants. The stress resultants are determined by substituting Eq. (3.44) into the 
following integral 


({JV}, {M}) = f 2 {o-} (1 ,z)dz (3.54) 

/2 k 

Substituting Eqn. (3.34) for the k?h layer stresses in the above equation and performing 
the necessary integration leads to the stress resultants of a composite laminate panel as 


'N\ = [M 

Mi L[s] LDllUJ K 


(3.55) 


where [A], [B], and [D are the extensional, coupling and bending stiffness matrices of 
the laminate, respectively, which for an n-layer laminate are defined as 

M=E[s] (z* + i~ z *) ( 3 - 56 ) 


I 
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(3.57) 


[5]=z|[e] (zL-4) 

k - 1 ~ k 

[D]=t}M{4„-^ (3-58) 

k = 1 * k 

The force and moment vectors resulting from the piezoelectric effect are defined as 

(RIK})= (tpW\ k E }l (\,z)dz (3.59) 

3.5.3 Stress Resultants for Small Deflections 

The piezoelectric force and moment vectors will be subsequently examined in much 
greater detail in Section 3.5.5. Nevertheless, the overall force resultant vector may be 
expressed in terms of the nodal d.o.f.’s, and are given here for the linear small 
displacement approximation as 

M = Me. ][r m ]{ w m } + [flic, Jr, ]K}- {tf, } 
=wsj^}+[£pj{n}-K} 

= R} + W-K} (3.60) 

Similarly the resultant moment vector may also be expressed in terms of the nodal d.o.f.'s 
as 

{M} = MC. R ]{w„ } + [D\C b \T b ]{ n } - {w, } 

=[sRK}+M8 ( W-H} 

= {Mj+{M b }-{M ,} (3.61) 
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3.5.4 Stress Resultants for Large Deflection 

The von Karman large deflection strain-displacement relations of Eq. (3.34) are 
substituted into Eq. (3.55) to determine the resultant force vector, hence 

M=Mk}+Mk°}+WW-K} 

= [^][cJrJK} + tMMcjr i ]K} 

+Mcjrj n }-K} 

= [ A I B „ ]K } + j l A \ e \ B o Kn } 

= {N m } + {N B } + {N t }-{N,} (3.62) 

Similarly the resultant moment vector may be determined as 

= [B\C m lT m ]R } + l - [EMC, \T„ ]R } 

+MQK]rHr} 

= + PPJR} - {M,} 

= {Mj + {M B } + {M i }-{M,} (3.63) 

Comparing Eqs. (3.62) and (3.63) with Eqs. (3.60) and (3.61), the resultant stress {jV^} 
and {M b } are the components due to large deflections. 
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3.5.5 Piezoelectric Resultant Forces and Moments 

Since each piezoelectric layer contributes to the total force resultant vector a 
summation over the range of np piezoelectric layers must be incorporated to account for 
each piezoelectric layer. Thus the piezoelectric force resultant vector is defined as 

W-ZI [Q]{d} k E n dz (3.64) 

k = 1 ** k 

The lamina stiffness [q], stress/charge constants {d}, and the electric field {E 3 } remain 
constant for each piezoelectric layer with respect to the 5-direction and are also assumed 
to be isotropic in the 1 - and 2- direction, thus the indicated integration reduces to 

k = 1 

= M{d\\ ••• \Ql\d\A - \q\M}«pK\e,'i (3.65) 


Furthermore Eq. (3.37) may be substituted for {E 3 } producing 


{N,} = -lQ\{d)A - [Q\{d} k h k ... {q\MUK\ b M} ( 3 - 66 ) 

K} = -K][sJK) (3.67) 

where 

- [Q\{d} t h t -. [Q], p {d} np h np \ (3.68) 

Similarly the piezoelectric moment resultant may be expressed as 


K}=-[^kK) 


(3.69) 


where 


[p u ]= 


^[Q\{d)A{z 2 + z,) ••• ^[Q\{d} k h k (z M +z k ) 


1 r^r 


,p{d}„ p Kp( z »^ +z „ P ) 


(3.70) 
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3.6 Element Matrices 


3.6.1 Introduction 

Evaluating the terms in the Hamilton's variational statement of Eq. (3.53) results in 
the finite element matrices. Variation of the potential and kinetic energies leads to the 
development of the element stiffness and mass matrices, respectively. During this 
investigation, body and concentrated forces are neglected. 

3.6.2 Linear Stiffness Matrices 

The finite element linear stiffness matrices will be determined first by evaluating the 
potential energy terms of Hamilton's variational statement Eq. (3.53). Thus the potential 
and electrical energy terms for the ifl 1 layer of Eq. (3.53) may be expressed as 

\({Ss} T {D})W 

¥■ 

= j[ft({&°} r +z{SK} r \a} k dz- J* (3.71) 

where the first integrand represents the strain energy and the second integrand represents 

the electrical energy due to the polarization properties of the piezoelectric in the 3- 
direction only. Thus by applying the stress resultants of Eq. (3.54), the variational energy 
becomes 

&(SE 3t )D 3k dz]dA (3.72) 

if- /2 J 

The first term of Eq. (3.72) is evaluated by substituting the force resultant vector of Eq. 
(3.60) which yields 

A A 

A 

+ [B\C b \T„ ]{w b } - to ]{E 3 })dA (3.73) 
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Thus by applying the results of Eq. (3.73), the first term of Eq. (3.72) may be expressed 
in the form 


\{9w m } T [B m f {[A\B m } + [Bp, K } + [P N % K ( 3 - 74 ) 

A 

Similarly by substituting the moment resultants of Eq. (3.61), the second term of Eq. 
(3.72) becomes 

\{Sk} T {M}dA = ([B]{s° } + [DJk) - {M, })dA 

A A 

= f [T b f [ c„ Y ([B][c„, \T m ]{w m } 

A 

+ [D\C b P K } - [P u ]{E U })dA (3.75) 

By collecting terms, Eq. (3.75) can be simplified as 

J{ } T [B b Y ([B\B m ]{w „, } + [Dp, ]{w „ } + [P u ][B, })dA (3.76) 

A 

The third term of Eq. (3.72) is evaluated by substituting Eq. (3.47) and Eq. (3.14) for the 
electric displacement density as 

J0(<® 3 jD 3 ,<fepl = 

A 

ft (3-77) 

A lk=\ k J 

Again integration of the piezoelectric lamina with respect to the thickness may be 
simplified by considering the geometric material assumptions previously mentioned. 
Thus Eq. (3.77) reduces to 

ft[(ffi* Ml [q\ K {®* }+ Ml \Ql + z, ) 

- (5E ik [ei {. d} k E lk h k + (SE 3k )e» , E u h k JjfiA (3.78) 


31 



By considering the laminate and utilizing matrix notation, Eq. (3.78) may be expressed as 


A 

+ {^} r K] r ([e»]-|/])K}]^ (3.79) 

where the matrices and [y\ are defined as 



WffelMi ”• 0 

M= i {d\llQ\{d} k J (3.81) 

0 - {d}l p [Ql P {d}„ p 

Thus combining Eqs. (3.74), (3.76), (3.79) yields Hamilton's variational statement, less 
the kinetic energy terms, which can be expressed as 

\\{Sw m } T [Bj[AlB m lw m } (3.82a) 

A 

+ (3.82b) 

+ {Sw„Y[B m \ T [PABM} ( 3 - 82c ) 

+ {Sw b Y[B b ] T [BlB m lM, m } (3.82d) 

♦ (3.82e) 

+ {5w i } T [B b ] T [P u \BM} ( 3 - 82f ) 
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(3.82h) 


+ {*»„ ) T [B,] T [P u ] T [B b ]{w b } 

+ K } r KNr]' - [e» Dh = ' 0 

The element stiffness matrices are determined from Eq. (3.82) and 
matrix form of 


'Sw b 

r 

[h\ [K,A foj 

V 

<K, 

> 

foj fo,] foj 

m 

/*V 


.M Ivl foL 

L <p ) 


where the corresponding linear stiffness matrices are 

[*.] = \[bJ [A\B m ]dA 

A 

[*«,]« \[Bj[B\B b ]dA 

A 

M= 

A 

[*.]= J[B.f 

M- Jtof 

A 

K]=\b*1[pA t \[KW 

A 

M=[sJ r (M-k 3 ]M 


(3.82i) 
may be cast into a 

(3.83) 

(3.84) 

(3.85) 

(3.86) 

(3.87) 

(3.88) 

(3.89) 

(3.90) 

(3.91) 

(3.92) 
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3.6.3 Large Deflection Equations 

The nonlinear element stiffness matrices are determined by following the same 
procedure outlined for the linear stiffness matrices in Section 3.6.2. In order to determine 
the nonlinear stiffness matrices, the von Karman large deflection strain must be included 
in the Hamilton's variational statement. The resulting variational potential and electrical 
energy statement, including the von Karman large deflection strain, may be expressed as 



np -r 

Ef 

*=1 k 


[{SE ik ){{d} T k [Ql{t 


}+z{/c} {d} k E ik + E ik 6 3ik |jrfz]c64 (3.93) 


Considering the stress resultants of Eq. (3.54) and substituting the results into the first 
integrand of Eq. (3.93) yields 


j[{&“ f {n}+ {<5*r} r {l M }}lA 


(3.94) 


A 

Since the membrane strain and curvatures are independent of the plate thickness, the 
variations in Eq. (3.94) becomes 


{&■ F = } r K Y + {<H } r [Be Y m T 

(3.95) 

{SkY ={Sw b } T [Bj 

(3.96) 


where the following relation is utilized 


4 \ [OlB, R }1 = ~ [SG\B e ]{w b } + i [e\B e }{5w b } 

J A l 


={e\B e lSw b } 


(3.97) 


The second integrand of Eq. (3.93) may be similarly evaluated by substituting the von 
Karman large displacements of Eq. (3.34) and performing the variational operation 
described above in Eqs. (3.95) and (3.96). Evaluation of the linear terms follows the 
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procedure outlined in Eqs. (3.78) and (3.79). Hence the variational statement of Eq. 
(3.93) becomes 


\[{Sw m [B m Y ]R } + \ [A\0\ B s ]K } 

+[s][sJK}-K}) 

+ {<K } r [B e Y [0f (Ms. ]K , } + \ Uim In } 

+ [s][sJK}-K}) 

+ {<K Y [5, r ([SR ]K } + \ [B\6\B e ]K } 

+MsJK}-K}) 

+ {*-, } r [B t Y fa Y (k K ) + \ Ms* IK }] 

+ {R} 7 [sJRF[sJK} 

+ Ik } r [s„ F tr] - [4 ]K }]^ = 0 (3.98) 

The variational statement contains identical terms which lead to the linear stiffness 
matrices, however nonlinear stiffness matrices resulting from the von Karman large 
deflection will appear and are indicated by the inclusion of the slope matrix. Thus the 
variational potential energy statement in Eq. (3.98) may be written as 
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Note that the linear expressions indicated in Eqs. (3.99a, c, d, i, k, m, o, p) are identical to 
those of Eqs. (3.82a, b, c, d, e, g, h, i) respectively. The remaining expressions which 
contain the slope matrix will lead to the nonlinear stiffness matrices. The following 
transformation relationship is applicable for any force resultant vector and will be utilized 
to further simplify expressions for the nonlinear stiffness matrices. Thus the product of 
the transpose of the slope matrix and any force vector may be expressed as 
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(3.101) 


Thus by applying Eqs. (3.9) and (3.32), the above relation yields 

[Of {N , } = [N, m = [N, \C 9 ]{«} = [N, \B e ]{w b } i = b, m, 0, <j> (3. 102) 

Equations (3.99e) and (3.99h) may be manipulated into a symmetrical form. Thus Eq. 
(3.99e) becomes 

\[BA T [e] T [A\B m ]{w m }dA = 

A 

k^y ji« # ntfrpp.k}+KRk}j« (3- io3 > 

^ A 

and Eq. (3.99h) becomes 

{Sw b Y \[B 6 i[e} T [P K \B,]{w,}dA = 

A 

ksw b Y Jfo f ([tff [P N K K }-K] [B, K }>« (3 • 1 04) 

** A 
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where the membrane force vectors of Eqs. (3.62) and (3.67) are used. Combining Eq. 
(3.103) with Eq. (3.99b) the resulting expression becomes 

[B e } r W{A\B m }{ Wm } 

A 

piKi 

+ \ {8"J [B a Y ]{w„ }^dA (3.1 05) 

The above relationship may be expressed in terms of the first order nonlinear stiffness 
matrices as 


| ({<H } : r ]{' w m } + {Sw b } T [n\ Nm ]{w s } + {dw m } T [n\ mb ]{■ w b }) (3.106) 

where the first-order nonlinear element stiffness matrices are given by 


KJ= \[B e M[A\B m W 

A 

(3.107) 

KJ = peUNABeW 

A 

(3.108) 

KJ= \[Bj[A\elB g }dA 

A 

(3.109) 

Equation (3.104) may be expressed as 

^ ({<H Y [«1( * ]K }+{dw b f \n\ st \w b }) 

(3.110) 

where 

KJ= \{B e i m 7 dA[P N ][B,} 

A 

(3.111) 

[«1 J=-\[B e i[N,][B e \U 

(3.112) 


A 


Similarly Eq. (3.99g) may be expressed symmetrically utilizing the bending force vector 
and adding it to Eq. (3.99j), hence 
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J \{S^} T {B e ] 7 m T {BlB b l Wb } 

+ f [5, f {B\9\B e fw b }}M (3.113) 

Equation (3.113) can be expressed as the first-order nonlinear stiffness matrix due to the 
laminate coupling matrix [5] and large deflection effects as 

^{Sw b } T [n\ m ]{w b } (3.114) 

where 

[n\ m ] = j \b, Y V>Y [B\B b ] + [B, Y [N t \B S ] + [B b Y [B\e\B $ ])« (3.115) 

A 

The first-order nonlinear stiffness matrix due to electromechanical coupling can be 
determined from Eq. (3.99n) as 

(3.H6) 

where 

[«1 (3-117) 

A 

The second-order nonlinear stiffness matrix is determined from Eq. (3.99f) as 

^{Sw b } T [nl b ]{w b } (3.118) 

where 

\[B e i[9YUl6\B e W (3.119) 

combining all the stiffness matrices the complete variation statement becomes 
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3.6.4 Mass Matrices 

The variation of the kinetic energy terms of Eq. (3.43) leads to the element mass 
matrices, which consist of mechanical quantities only, since the electrical d.o.f.'s do not 
have an equivalent inertial analogy. Thus the variational kinetic energy may be written as 


J p{^dw} T {w} + {du} T {u\ + {Sv} T {v})dV 

v- 

- - 1 p{^dw} T {w} + {Su} T {ii} + {8 v} t {v})c/F 
v- 

= 4{^} r KK^} + {* 4, m } r K]{w„}) (3.122) 

where the element mass matrices are given by 

KHZir \{H w }hp{x,y\_H w \dA[T b \ (3.123) 


[m„] = [T m ] T \ \{Hjhp{x,y\HA+\{H v }hp{x,y\HA]dA[T m } (3.124) 


3.6.5 External Force Vector 

In completing the variations indicated in the Hamilton's variational statement of Eq. 
(3.43), the work done due to external forces, body forces and surface traction's were 
assumed negligible, however the electrostatic work due to the externally applied voltage 
must be included. The electrostatic work done as, described in Appendix C, is given by 
Eq. (C.23). Since the space charge within the volume of the piezoelectric is zero, only a 
surface charge accumulates on the electrode surfaces, hence the virtual electrostatic work 
done can be expressed as 

\3V Pr JS = J \ 8 \V 4 , } r {p cs }dA = -{dwp } r {p ^ } (3.125) 

S 2 S 2 


where 

\p*} = -\{PoW (3.126) 

A 
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3.6.6 Element Equations of Motion 


The element equations of motion may be formed by substituting Eqs. (3.120), (3.121), 


and (3.125) into Hamilton's variational statement Eq. (3.53), and collecting terms. 
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(3.127) 


3.7 System Matrices 

The element equations of motion Eq. (3.127), are a set of equations which describe 
the fully coupled structural and electrical properties. Application of Eq. (3.127) requires 
the implementation of an assembly procedure in accordance to the prescribed electrical 
and structural boundary conditions. The assembly process for the structural stiffness can 
be shown symbolically as 


[*] = £[*] (3-128) 

where the global stiffness matrix [AT] has dimensions m x m for m structural and electrical 
d.o.f.'s and the element stiffness [&] is of size (24+np)x(24+np). The assembly procedure 
can be visualized by first starting with a null global stiffness matrix, then subsequently 
adding to it [&] of each element until ail the elements are considered. Assembly of the 
mass matrix is accomplished using an identical procedure, however special attention is 
required for the piezoelectric elements and will be subsequently discussed in greater 
detail. Assembly of Eq. (3.127) yields the following fully coupled system of equations 
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(3.129) 
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If the given structure contains several piezoceramic patches, and each patch consists 
of n finite elements (Figure 3.4), the assembly process must be modified only for 
elements which contain piezoelectric material. The prescribed electric boundary 
conditions require that the electrode be maintained to an equipotential, therefore each 
patch must consist of one electrical d.o.f. and can be simply assembled as 

KMk), - kl- kU r ( 3 - 14 °) 

where N is the number of patches and {w^ } as defined in Eq. (3.35) is a npxl vector 

representing np number of piezoceramic layers. The solution for the assembled system of 
equations in Eq. (3.129) may be obtained by utilizing the standard finite element solution 
procedure. Thus during the solution process there is no need to distinguish between 
structural or electrical quantities other than known or unknown quantities. 

3.8 Solution of Static Sensor Equation 

To determine the voltage produced by a single piezoelectric patch bonded to a panel 
subjected to a static uniform distributed load, Eq. (3.129) may be partitioned as 




[[Kj+\[m^]j{W} + K,W,= 0 
where Eq. (3. 142) may be solved in terms of the unknown voltage as 


(3.141) 

(3.142) 


W^-K;'{[K^]+±[N\Jyw} (3.143) 

where {P ^ } = 0 since there is no externally applied voltage to sensor. Furthermore Eq. 
(3.143) may be substituted into Eq. (3.141) resulting in a system of equations which may 
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be solved for the structural d.o.f. in response to the structural loading. The structural 
d.o.f. may be subsequently applied to Eq. (3.143) in order to determine the voltage. Since 
the system under consideration is static, the inertial terms will be identically zero. In 

addition since there is no externally applied voltage, the first-order nonlinear electrically 
coupled stiffness [TVIjyJ will also be identically zero. Substituting Eq. (3.143) into Eq. 

(3.141) yields the following system of equations, which must be solved for the structural 
d.o.f. due to the static loading 

J-^^KIKIK 0 '} = iKi ( 3 - 144 ) 

The Newton-Raphson iterative method is used to solve the nonlinear system of equations 
of Eq. (3.144). Recall that the first- and second-order nonlinear stiffness matrices are 
functions of the unknown displacements. 

3.8.1 Newton-Raphson Iteration Method 

The Newton-Raphson [20] iterative method is used to solve Eq. (3.144) by solving for 
an incremental deflection which is given by 

m M ={W} l + {AW}, (3.145) 

The iterative procedure is carried out until the incremental deflection {A W} approaches 
zero resulting in a converged static deflection {W } . 

Equation (3.144) can be written as a function, which may be expressed in terms of a 
Taylor series expansion thus 
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(3.146) 


w}=a«j+j [ ^2]+i [ M] 
-KAK,Y[K^[K H lK,]'[mJ 

The truncated Taylor series expansion is represented symbolically as 

{' ¥(W + AW)} = {^(IF)} + | ^^H {Ar} = {0} (3.147) 

Differentiation of Eq. (3.146), referred to as the tangent stiffness, results in 

cW(W) = 4(K] + ^2] + i[M] 

- K, K l 1 M - ~ K, R 1’ i m J 

-\[NUW KT K l 1 K (3-148) 

The indicated differentiation of the tangent stiffness in Eq. (3.148) is accomplished by a 
term by term evaluation, beginning with the first term, which is the linear structural 
stiffness matrix. Since the linear stiffness matrix is constant with respect to the structural 
d.o.f., the differentiation is a trivial operation resulting in 

d([KjW}) = [Kj{dW} (3.149) 

Differentiation of the first-order nonlinear stiffness matrices invokes a two-step approach 
since [iVl] is comprised of and [M ra ]. The first step will deal 

with [AT*,,, ], [M A/ji ], and [AT tol ] , and the second step will involve the differentiation of 
[N1 m ]. The differentiation must be performed on the element level stiffness matrices, 
thus [n\] is expressed as 
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(3.150) 


(M)M 


fa*,.] t «i&„r 

IKJ [o]. 


w. 


w„ 


Thus the derivative of [Nl] can be performed on the partitioned element level definitions 
given in Eq. (3.150) as 

d (R,„ ]{ }) + d (K„, ]{ w „, }) (3.151) 

and 

4RJH!) (3.152) 

Applying the definitions expressed in Eqs. (3.107) and (3.108) to Eqs. (3.151)and (3.152) 
yields 

”(RJ+ J[c,„fUPlc e ]{oR 

+ Pir J[Qf (3.153) 

where [Q ]{«} = {6 } . Performing the required differentiation of Eq. (3.153) yields 


\[cJ[Al{dem+[elde})dA 
+[Tj \[C 6 i[ddUA[C,MdA 

A 

+ [Tj \[CeiM[A\C m ]{db}dA (3.154) 

A 

Simplification may be achieved in the above equation by considering the following 
relations 


\do]{e} = 


dw, x 
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dw, x w, x 


dw, y w„ 

dw, y w, x +dw, x w, 


(3.155) 


y J 
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°lr rfw i f dw ” w ” 

[0]{d6) = 0 w, y \ dw X \ = < dw >y > (3.156) 

_w, y w , J >yJ [dw, y w, x +dw, x w, y ^ 

Therefore 

[dd\o} = [e]{de) = [elc e lda} ( 3 . 157 ) 

Thus Eq. (3.154) may be cast in the following form 

[Tj \[cJ[A\elC a \iA{da} 

A 

+ {Trf \{C e i[N m lC e lda}dA 

A 

*[T b J \[C s i[e] T [A\C m W{db\ (3.158) 

A 

Thus the differentiation procedure can be summarized as 

] + K }| + d{^[n\ mi ]{w b } ^ 

= l dw b } + [«V, K<M } + [n\ bm \dw m } (3.159) 

The second step used to differentiate [NI nb ] uses a similar approach involving Eqs. 
(3.114), (3.99g), and (3.99j) resulting in 

= [Tj \[Ceim T [B\C b Ib}dA 

^ A 

+ \[Tj J[C t f [BMcMdA (3.160) 

Performing the indicated differentiation yields 

[Tj \[C e i{deT[B\CMdA + [Tj \[C e f [0] T [B][C b ]{da}dA 

A A 

+ [T b ] T \[C b i[Ble]{de}dA (3.161) 
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and by using the relations in Eq. (3.157) yields 


[Tji([c e f[N b ic e ] + [cjmm] 

A 

+ [C b f[B\elC e ])dA{da} (3.162) 


which may be summarized as 



[nl m ]{dw b } 


Combining the results of Eqs. (3.159) and (3.163) yields 


d\ h/iajir}) = [N\]{dw\ 
2 J 


(3.163) 


(3.164) 


Similarly, differentiation of the second-order nonlinear stiffness matrix may be performed 
using the element expression in Eq. (3.119) and substituting [C 0 ]{a} = {#}, thus 

\[n2 b ]{w b } = hrj \[C e i[d] T [A\0]{e}dA (3.165) 

Performing the indicated differentiation operation produces 

= far jtaf mmKWM = o.i66) 

which may be expressed as 


( 

d 

\ 


±[N2]{W}^ = lN2]{dW} 


(3.167) 


The next term in Eq. (3.148) involves differentiation of the linear electrical-structural 
coupling stiffness matrices, which is a trivial operation since the stiffness matrices are 
constant with respect to the structural d.o.f., hence 
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= KA K X\ K J. dw ) 


(3.168) 


The next term of Eq. (3.148) involves differentiation of the nonlinear electrical-structural 
coupling stiffness matrix, which is a function of the bending displacement, thus the 
differentiation is indicated as 




(3.169) 


which can be expressed in terms of the element level matrices as 


TJ' 

V^lMj 


hr 


Kiry 

L to] JJi 


(3.170) 


A w .. 


by recalling Eq. (3.117), the first order nonlinear coupling stiffness term is defined as 

= ftfftr jMQMrJKJ (3.171) 

^ A 

since it was previously shown that the linear stiffness matrices are constant matrices they 
are not affected by the differentiation indicated in Eq. (3.169). Thus differentiation may 
be accomplished by differentiating Eq. (3.171), which results in 

[b„ r ft r \[M\c e wiT,, ]{w t } + ft r ft r jmq mt> j o . i 72) 

A A 

by using the fact that {a} = \T b \w b }. The relations of Eq. (3.157), Eq. (3.172) may be 
reduced to 

[sJftrjMQMrj (3.173) 

A 

Thus the differentiation indicated in Eq. (3.169) can be summarized as 

KMtV WJ.dW) (3.174) 
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The next term of Eq. (3.148) requiring differentiation is the matrix transpose equation of 
Eq. (3.169) 




(3.175) 


thus neglecting the linear stiffness matrices differentiation of the first order nonlinear 
coupling matrix is simply the matrix transpose of Eq. (3.172) and is given by 

[Tj J[C„f [de\ r dA[P N iBM) + jpT jtaf [efdA[P N }[B,] (3.176) 

A ^ A 

thus differentiation results in 


= [ N Kj^\'[KjdW} (3.177) 

The next term of Eq. (3.148) involves the product of two nonlinear coupling matrices 


d[^[N\ w ,\K,]'[N\Jw} 


(3.178) 


which is expressed at the element level as 
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k«Jl 


\KV 


KiJlfn 
L [0] Jk 


(3.179) 


4L [0] J 

This results in a differential with respect to the bending displacement only, resulting in 


d {^ rA ti\ k X Ikk}] (3.180) 

where 

kJ=Kf \[C e im T dA[PAB,]= Kr (3-181) 

A 

Once again the linear stiffness matrix is constant and unaffected by differentiation. Thus 
the differentiation is performed on the product of the first-order nonlinear stiffness 
matrices only resulting in 
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UtJ j[c 9 f[deTdAlP„tB,h\%Y{p N J \[e\c e ]dA[TM;\ 

4 A A 

+\[T b ] T \[c e i[9\u[p N \B JkJ%Y[p N Y \[e\c e Mr b ]{w b } 

4 a A 

+ j[tJ Jfc f [OW[P, \k t ]% ] [P N F \[8\C 6 \1A[TA (3.182) 

4 A A 

which is simplified by employing the relations of Eq. (3.157) resulting in 

[n\ H ] [k, ]-' [nl „ K }) = | K ] h ]"' KK) (3-183) 

Reassembling the system matrices yields 

d{^[NKA K *Y[ m t»\w}) = \[m w ,\K,] l [mJ_dw} 0.i84> 

and collecting all terms yields the tangent stiffness matrix given by 

= [KJ + [N2}+ [m]-[K^\K,Y[K^] 

- K* K ]'' fa*. ]- [M J[K, P K 1-f i m * l K * f KJ (3.185) 

Thus Eq. (3.147) becomes the nonlinear static equation where 

[K,J{AW} = {APJ (3.186) 
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The incremental right-hand side vector is 


{A/D = tf,}-(K]+ j[JV2]+i[M] 

- K,W[k^}- 1 -{k^\k, f'KJ 

“1*1. (3.187) 

In the Newton-Raphson iteration, the initial nonlinear stiffness matrices are determined 
from the linear static displacements. Once the incremental load vector and tangent 
stiffness matrix are assembled, the Newton-Raphson tangent equation shown in Eq. 
(3.185) is solved to determine {AP W } , which is subsequently used to update the static 

deflections. During the iterative process, the incremental load vector and the incremental 
deflection will approach zero, resulting in a converged deflection solution. 
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Fig. 3.3 Isotropic panel with surface mounted piezoceramics 
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Fig. 3.4 Piezoceramic elements for two patches 
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CHAPTER IV 
EXPERIMENTAL SETUP 
4.1 Introduction 

It was desired that a typical representative aircraft structural panel be tested, and the 
experimental results compared to the finite element analysis. Therefore a 14 n xl0 M x0.040" 
aluminum panel was selected for the experiments. For a panel of this size, it was much 
easier to construct a test fixture that achieves approximate clamped boundary conditions 
as opposed to simply supported boundary conditions. 

4.2 Establishment of Clamped Boundary Conditions 
A clamping fixture consisting of two aluminum plates of size 25"x20"xl", with a 
rectangular 14"xl0'’ hole machined out of the center, was manufactured to achieve the 
clamped boundary conditions as shown in Figure 4.1. The rectangular cut-out was 
centered such that there was a minimum of four inches of clamping surface around the 
perimeter of the 14"xl0" panel cut-out. The mating surfaces of the clamping fixture were 
measured to be within ±0.002” flatness. A series of 26 pilot holes for 1/4” cap screws 
were drilled in an evenly spaced pattern 1/2” beyond the perimeter of the panel cut-out. 
Another series of 8 pilot holes were spaced evenly around the outer perimeter of the 
fixture 3/4” from the edges and sized to accommodate the 3/8” cap screws. The 0.040" 
test panel specimen was cut to 25"x20” so as to fit the overall clamping fixture. Two 1/4” 
alignment pins were placed at opposite comers to assure repeatable and consistent 
assembly. Wooden support blocks were placed on the bottom of the panel to support the 
clamping fixture in a vertical fashion. The clamping fixture was assembled with grade-8, 
high strength hex head cap screws, washers, and nuts. The 1/4” and 3/8” cap screws were 
tightened to 15 and 20 ft-lb respectively in 5 ft-lb increments. The tightening sequence 
followed is depicted in Figure 4.2. 
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4.3 Piezoelectric Wafer Preparation 

A single 2.5”xl.5"x0.010” piezoelectric wafer was bonded to the panel 3/4” from the 
boundary, symmetrically on the centerline (see Figure 4.3). Before the wafer can be 
attached to the panel, electrical leads must be attached to the wafer electrodes to facilitate 
the instrumentation. Since the piezoelectric wafer was polarized in the 3-direction the 
electrodes are on the top and bottom surfaces. The piezoelectric is a polarized dielectric, 
therefore if any residual electrode material is attached to any of the edges in the 2- or 3- 
direction, a complete or near short circuit will adversely affect the piezoelectric process, 
if not completely inhibit it. In addition, since the electric field is the gradient of the 
electrical potential as shown in Eq. (3.36) a very large electric field, or potential gradient 
will exist. This condition could cause dielectric breakdown of the piezoelectric resulting 
in a short circuit condition. This adverse condition may be minimized if the piezoelectric 
edges are straight and perpendicular. To assess the quality of the piezoelectric edges the 
resistance and capacitance can be measured with a high quality multimeter and 
comparing the results to the manufacturer's specifications. If an ill-conditioned edge is 
detected, as indicated by a low resistance or capacitance, then the edges may be gently 
scraped or sanded with a very fine machinist wet abrasive paper. If the specified 
resistance cannot be obtained then the piezoelectric wafer should not be used. 

4.4 Piezoelectric Lead Attachment 

Electrical leads must be soldered to the electrodes of the piezoelectric wafer. 
Traditionally to obtain a good electrically conductive and mechanically reliable solder 
bond both the substrate and the lead being attached must be raised to a temperature 
greater than the melting point of the solder. This technique will adversely affect the 
piezoelectric wafer since depolarization will occur if subjected to temperatures in excess 
of the Curie point which is 365°C for PZT-5A [21]. The depolarization may occur 
locally within the piezoelectric, thereby altering the overall charge constants. A complete 
depolarization rendering the piezoelectric inoperative is unlikely, therefore extreme care 
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must be taken to assure that the piezoelectric properties are not altered during the 
soldering process. 

The piezoelectric wafer is a polarized dielectric, and since thermal conductivity is 
proportional to electrical conductivity [22], the piezoelectric wafer is a thermal insulator. 
Hence the heat applied during soldering will not readily dissipate through the wafer. On 
the other hand the electrodes are nickel and are highly conductive, however since they are 
only 0.00005 - 0.0002” thick their thermal capacity is very small. The only other thermal 
member remaining is the lead to be attached. The lead size should be selected to 
minimize its thermal capacity so that it does not act as a heat sink and direct the heat 
away from the wafer, thus increasing the soldering time. The leads used were copper foil 
2"xl/16' , x0.00r\ A temperature selectable thermostaticly controlled soldering iron fitted 
with a small pencil tip was used and set to 360°F. The electrodes of the piezoelectric 
wafer were cleaned with alcohol and a small amount of liquid flux applied. The lead was 
cleaned, fluxed, and tinned with solder. After the lead cooled additional flux was applied 
to the lead and placed on the substrate, then with light pressure, a tinned soldering iron 
was placed on top of the lead until the solder between the lead and electrode melted. The 
soldering iron was removed within five seconds as recommended by the piezoceramic 
manufacturer. 

4.5 Piezoelectric Wafer Bonding 

The piezoelectric wafer attachment is carried out in a two step process to ensure both 
electrical isolation and mechanical bonding. A strain gage epoxy adhesive system was 
used to bond the wafer to the aluminum panel [23]. The electrical isolation is obtained by 
applying a thin layer of adhesive to the panel prior to the application of the wafer. 
Initially the panel is cleaned with alcohol to remove any dirt and oil. The piezoelectric 
location was measured and marked within 1/64” using a sharpened mechanical lead 
holder. The scribe lines were extended three inches beyond the actual location for future 
reference since excess adhesive will cover the lines as indicated in Figure 4.4. Mylar 
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adhesive tape was applied 1/16" beyond the scribe lines creating a rectangular mold 
which will be filled with adhesive. Additional mylar tape was used to cover and protect 
the scribe lines remaining outside the frame, see Figure 4.5. The framed area was wet 
sanded with alcohol and abrasive cloth, wiped clean and treated with a micro-abrasive 
metal conditioner. Next a semipermiable Teflon cloth and mylar sheet were cut and 
placed over the area to be glued and one edge was attached to the panel with mylar tape. 
The adhesive is mixed and a thin layer is applied to the framed area, the Teflon cloth was 
placed over the adhesive and covered with the mylar sheet and a vacuum pad was 
attached and activated. The vacuum pad will provide a uniform pressure permitting 
excess adhesive to penetrate the Teflon cloth. The adhesive was permitted to cure 
overnight since the minimum curing time was six hours. 

Once the adhesive cured the vacuum pad was removed along with the mylar tape. 
The Teflon cloth provided an abraded surface texture however the adhesive should be 
uniform without any voids or bumps. Using the original exact scribe lines reapply the 
mylar tape to create a frame which will be exactly the size of the piezoelectric wafer. The 
adhesive surface and the piezoelectric wafer were cleaned with alcohol and a thin layer of 
adhesive was applied to the existing base coat on the panel. The piezoelectric wafer was 
placed on the adhesive and covered with a thin layer of foam and the vacuum pad placed 
over it and allowed to cure over night. It should be noted that the electrical lead on the 
wafer should be as thin as possible including any residual solder. If however the 
piezoelectric wafer is very thin, and subjected to an excessive vacuum loading, the lead 
may crack the piezoelectric wafer. For a flat horizontal application, only a slight vacuum 
loading is required, to ensure a uniform distributed load during the curing process. 

4.6 Uniform Distributed Loading 

To assess the static piezoelectric sensor a uniform distributed load was applied to the 
panel. However, since piezoelectrics cannot sustain a static charge the uniform 
distributed load must be applied instantaneously. To achieve the instantaneous load the 


61 



test clamping fixture was modified to include a 14"xl6"xl/4" plate on one side, thus 
creating a vacuum chamber. The vacuum plate was attached using the existing fixture 
bolts and was sealed with a vacuum grease. Three fittings facilitated a vacuum pump, 
vacuum gage, and a quick-release ball valve as indicated in Figure 4.6. To conduct the 
static sensor experiment the vacuum was applied and the maximum plate deflection was 
measured using a dial indicator. Once the desired deflection was achieved the vacuum 
was held by closing the valve and noting the corresponding pressure. With the plate in a 
deformed state, the residual charge was allowed to dissipate. At this time, the ball valve 
was opened quickly and the sensor voltage was measured by recording the time history 
using a fast Fourier analyzer. 

4.7 Material Properties 

The experiments conducted for this research were performed on a rectangular 
isotropic aluminum panel with a single bonded piezoceramic patch. Various static and 
dynamic experiments were conducted in order to obtain data that was compared to finite 
element analysis which will be discussed subsequently in Chapter IV. The piezoceramic 
used herein was PZT-5A manufactured by Morgan Matroc [21]. Table 4.1 provides the 
physical material properties. The physical properties of the aluminum panel used for the 
experiments may be found in Table 4.2. 
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PZT-5A 


Charge Constants 

d 31 

-171e-12 

m/V 

Permittivity 

€ r 

1700 


Density 

P 

7700 

kg/m 3 

Length 

L 

0.0762 

m 

Width 

W 

0.0381 

m 

Thickness 

h 

254e-6 

m 

Young's Modulus 

E 

6.9el0 

N/m? 

Poisson's Ratio 

V 

0.31 



Table 4.1 Piezoceramic Properties 


Aluminum Panel 

Density 

P 

2702 

kg/m 3 

Length 

L 

0.356 

m 

Width 

W 

0.254 

m 

Thickness 

h 

1.02e-3 

m 

Young's Modulus 

E 

6.1el0 

N/m 2 

Poisson's Ratio 

V 

0.31 



Table 4.2 Aluminum Panel Properties 
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Fig. 4.2 Bolt tightening sequence 


Fig. 4.3 Piezoceramic location 
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Fig. 4.4 Piezoceramic mounting preparation 
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Piezoceramic 


Fig. 4.5 Piezoceramic adhesive preparation 
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CHAPTER V 

NUMERICAL AND EXPERIMENTAL RESULTS 
5.1 Introduction 

The verification of the linear finite element model was conducted by considering a 
simply supported square panel. Initially, an isotropic panel is considered where the 
maximum transverse deflection was compared to the classical solution. Next a bonded 
piezoelectric patch is included and the predicted sensor voltage resulting from a uniform 
distributed load is compared to a voltage determined using a classical solution method. 
The large deflection finite element model was verified by comparing results to a single 
mode classical solution for moderately large deflections of a square clamped plate. 
Subsequently, finite element results are presented and compared to experimental results 
obtained from a clamped rectangular panel. Results include static and dynamic sensor 
applications and dynamic actuation. 

5.1.1 Static Small Deflection 

In order to verify the finite element model formulation, a simply supported isotropic 
panel subjected to a uniformly distributed load is considered. The small deflection 
approximation is used in the finite element formulation. The Navier solution for a simply 
supported rectangular plate is computed and used to validate the finite element model. 
The Navier solution for the transverse deflection of a simply supported plate under a 
uniform load is given by [24] 


'*(x,y)=- r z2,L- 

71 V m= \ „=} 


1 


mn 


f m 2 n 
ya 2 b 2 


2\ 2 

) 


■sin 


mrcc 


\ a 


sin - 


MDC 


V b 


(5.1) 


where m = 1,3,5, • • • oo and n - 1,3,5, • • • oo . Thus the maximum deflection can be found by 
substituting x = ^ and T = Figure 5.1 provides the finite element analysis non- 


dimensional maximum deflection compared to the classic Navier solution of a 
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10"xl0"x0.040" simply supported plate subjected to uniform loads. The accuracy of the 
finite-element, small-deflection formulation has thus been established. 

5.1.2 Static Large Deflection 

To verify the large deflection finite element model formulation a clamped square 
panel subjected to a uniform distributed load is considered. A single mode classical 
solution for moderately large transverse deflections is given by Chia [25] as 

/ \ 3 4 

+ 0.2522-^- = 0.0001333 isf— (5.2) 

h ) h Dh 

where w max =2.5223 w n , h is the panel thickness, D is the bending rigidity, q 0 is the 
uniform distributed loading, and a is the width of the panel. Comparison of the large 
deflection finite element formulation non-dimensional displacement and the single mode 
solution is shown in Figure 5.2. Excellent agreement was obtained for the large defection 
formulation, thus establishing the accuracy of the finite element model. 

5.1.3 Static Sensor 

In order to validate the static sensor formulation an isotropic plate with a single 
piezoelectric patch bonded to the surface is considered using a classical voltage solution 
method. The classic solution applies Gauss's law to determine the charge enclosed within 
a surface and is given as 


9 =jb<fe (5.3) 

where D is the electric displacement density and da is the differential normal vector of 
the surface S. However, the space charge within a dielectric is identically zero, thus Eq. 
(5.3) cannot be directly applied to piezoceramics. Since the charge generated by the k?h 
piezoceramic layer is accumulated at the electrodes (top and bottom surfaces of the kfh 
piezoceramic layer), Eq. (5.3) may be expanded to determine the effective charge such 
that 
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(5.4) 


<lk = 


1 


Jj" D 3 dxdy + JJ D 3 dxdy 


~\S l (z=z k ^) S 2 (z=z k ) J 

Substituting Eq. (2.12) into Eq. (5.4), setting the applied electric field to zero for the 
sensor application, and considering an effective or average electrode surface area yields 

q k = JjM(fc} + z{fc})dxdy (5.5) 

S' 

Substituting the strain-displacement relationships of Eqs. (3.16) and (3.17) yields 

q k = Jj(e 3 i(".,- zw .») + e 32 ( v k y -zw,„)+ e x {u, y +v, x -z2w, xy ))ixdy (5.6) 

Neglecting the inplane strain, since the small displacement approximation was used, Eq. 
(5.6) may be simplified as 

<lk = ff- 4 e il W ’*' +e 32 W ’v + 2 e 3 <i W ’ V } ixd .y ( 5 - 7 ) 

s e 

Utilizing the Navier solution of Eq. (5.1) for the transverse deflection and setting e 36 = 0, 
the charge of Eq. (5.7) becomes 


I6q g b ^ ^ v 1 
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(5.8) 


y y ) 

where Zp represents the distance of the mid plane of the tfh piezoceramic layer and is the 
effective electrode area as shown in Figure 5.3. The effective electrode area establishes 
the electrical boundary conditions to facilitate the electric displacement density. The 
classic and finite element voltages due to various pressure loading applied to a simply 
supported 10”xl0" plate is shown in Figure 5.4. The excellent agreement indicates that 
the finite element formulation will accurately predict piezoelectric sensor voltages. 
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5.2 Experimental and Analytical Comparison 

The experiment was conducted using a rectangular aluminum plate with a single 
piezoceramic patch bonded to the top surface. The validated finite element model was 
modified to include clamped boundary conditions. A full analytical model was used to 
model the complete experimental plate by using a 10x8 mesh which includes 4 
piezoceramic elements and is shown in Figure 5.5. By employing a high input 
impedance charge amplifier, the charge signal may be recovered and related to the 
voltage generated by the piezoceramic through the intrinsic piezoelectric capacitance as 
V=q/C , as described in Appendix C. Since a charge amplifier was used during the 
experimental test (refer to Appendix E) the charge was converted to volts during the 
analysis. It should be mentioned that if a rate of strain is desired, then a current amplifier 
must be used in lieu of a charge amplifier. 

The finite element simulation requires modal damping values of the plate, so they 
were determined experimentally along with the natural frequencies and are shown in 
Table 5.1. In the following Sections, a complete comparison of experimental and 
analytical results of the static and dynamic sensor and dynamic actuation is provided. 
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Experiment 

Analysis 

Frequency Hz. 

Modal Damping 

Frequency Hz. 

103.06 

0.034 

106.71 

164.56 

0.0047 

173.0 

236.1 

0.0095 

256.47 

268.81 

0.0007 

282.58 

298.94 

0.004 

316.70 

392.9 

0.0031 

422.65 

412.2 

0.0057 

433.65 

488.5 

0.0034 

491.73 


Table 5.1 Natural Frequencies and Damping Values 


5.2.1 Static Sensor 

The static sensor analysis was conducted in two steps. The first step implements the 
small deflection assumption, however in order to accurately predict the experimental 
results, a second approach including the large deflection approximation is performed. 
The linear piezoelectric theory accurately predicts the sensor voltage for small pressure 
loading, however for higher pressures the geometrical nonlinearities dictate the plate 
deflection and thereby affect the sensor voltage. Given the fully coupled electrical- 
structural formulation, the large deflection assumption includes the nonlinear electrical- 
structural stiffness matrices of Eq. (3.143). It should be noted however, that there is no 
nonlinear stiffness associated with the uncoupled electrical d.o.f. since the linear 
piezoelectric theory assumption is maintained. 

The small deflection analysis and experimental sensor voltages of a 14"xl0"x0.040" 
clamped plate due to uniform distributed loads are shown in Figure 5.6. Note that for 
small pressure loading, the small deflection analysis results coincide with the test results. 
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Thus for small pressure loading the small deflection analysis is valid. However it should 
be noted that for high speed flight vehicles it is not uncommon to experience differential 
pressures in the range of 3-5 psi, thus the experimental loading considered is not 
unreasonably high. The non-dimensional deflection due to the large deflection analysis is 
shown in conjunction with the experimental results of a 14"xl0"x0.040" clamped plate in 
Figure 5.7, and excellent agreement between the large deflection analysis and the test data 
is achieved. Similarly the experimental and analytical sensor voltage based on the large 
deflection analysis is shown in Figure 5.8. The agreement for the sensor voltage between 
the analysis and the measured experimental results is not as close as the maximum 
deflection. Deflections and sensor voltages from the small deflection analysis are clearly 
not agreeable with the test data. The static sensor analysis was based on the assumption 
that the uniform distributed loading may be modeled as a step function triggered at some 
time t 0 . Analytically, the piezoceramic is capable of producing a static charge, however 
real piezoceramics are not physically able to sustain a charge to a true DC response [26]. 
The static sensor experiment used was constructed to approximate a step response by 
establishing a distributed load and then allowing the generated charge to dissipate, then a 
ball valve was opened quickly thereby releasing the distributed load. Thus the pressure 
loading is actually a transient response with a finite rise time. As the distributed loading 
was increased the slope of the pressure discharge became more critical indicating that the 
discharge time could no longer be assumed to be instantaneous. Therefore improved 
analytical results may be obtained by simulating the transient response numerically and 
computing the subsequent response through Duhamefs integral. 

5.2.2 Dynamic Sensor 

The dynamic sensor analysis is based on the small deflection assumption since the 
applied loading was within the linear analysis range. However modal damping values are 
required and thus determined experimentally and are shown in Table 5.1. The 
experimental tests were conducted using a point load provided by an electrodynamic 
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shaker attached to the center of the panel as shown in Figure 5.9. A sinusoidal excitation 
signal was applied to the shaker and the subsequent piezoceramic voltage was measured 
and recorded. For the steady state dynamic sensor analysis the structural d.o.f. are 
transformed into modal coordinates by 

{W } = ^q r (t) {v)r = H{?} (5-9) 

r~\ 

where the q's are the modal coordinates and [y/\ is a reduced set of mode shapes. The 
small deflection equations of motion may be obtained by substituting Eq. (3.143) into Eq. 
(3.141) resulting in 

[M]{w}+(lKj- [e w Ie, [E„ Dpr } = {P w } (5.1 0) 

The equations of motion in Eq. (5.10) may be reduced to a set of uncoupled modal 
equations of much smaller d.o.f., by utilizing the modal transformation of Eq. (5.9), thus 
resulting in 

m r q r +c r q r +k r q r = f r r =1,2,...,« (5.11) 

where the small deflection assumption was imposed resulting in the linear equations of 
motion. The uncoupled equation of motion of Eq. (5.1 1) may be written in the following 
form as 

9, + 2 +a 2 r q r =— (5.12) 

m r 

The modal mass, stiffness and force are obtained from 

yf {W\[K-})W\ = (m r ,k r ) (5.13) 

where 

[K] = (K]- [k„Jk, ]-’ [k^ ]) (5. 14a) 

{r} T {Pj = fr (5.14b) 
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The modal damping £ ratios were determined experimentally, and are shown in Table 
5.1. Thus the modal coordinates qf in Eq. (5.12) may be determined and substituted into 
Eq. (5.9) to determine the nodal d.o.f., which are then substituted into Eq. (3.143) to 
determine the sensor voltage. Note that the first order nonlinear coupling stiffness of Eq. 
(3.143) is identically zero since the small deflection assumption has been applied. The 
dynamic analytical sensor voltage due to a sinusoidal point load applied to a 
14"xl0"x0.040" clamped plate is compared with the dynamic experimental sensor voltage 
in Figure 5.10. Similarly, the predicted and measured displacements of the plate center 
are compared in Figure 5.11. Excellent agreement between the experimental results and 
the analytical predictions are obtained for both the sensor voltage and the plate 
displacement as expected, given the linear static sensor results for small deflections. The 
dynamic analysis was computed using the reduced modal data set for the first eight 
modes shown in Table 5.1. 

5.2.3 Dynamic Actuator 

Since a single piezoelectric patch was used for the experiments, and it was observed 
that the resulting plate deflections were small due to piezoelectric actuation, the small 
deflection analysis was employed in the dynamic actuator formulation. In the 
experiment, a 120 Hz sinusoidal excitation signal voltage was amplified to 20 volts and 
supplied to the piezoceramic patch. The resulting acceleration at the panel center was 
measured as a time history on a fast fourier analyzer. The time domain acceleration 
signal was transformed in the frequency domain and integrated twice to determine the 
center plate displacement. The frequency domain integration was performed during post 
processing using Matlab©- 

The actuator equation can be determined from Eq. (3.129) as 

[M]{w} + (V J - i [N1 ]-[*„ I K t ]-' \kJ^{W} = -K» \K 0 ]-’ {/>, } (5.15) 
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where {P^} is the electrical loading determined from Eq. (3.126). Similarly the modal 
coordinate transformation shown in Section 5.2.2 can be applied to Eq. (5.15) resulting in 


™ r <lr+c r q r +k r q r =f r 

, f 

q r +2C r o) r q r +co r q r = — 

m r 

The modal mass, stiffness and force are obtained from 


(5.16) 

(5.17) 


={m r ,k r ) 


where 


K, \k, ]-' [k„ ] 


(5.18) 

(5.19a) 


-WY r [Kj[K,\'{P,\ = f r (5.19b) 

and f r is the applied electrical load determined in Eq. (3.126). Note that even though the 
small deflection assumption is employed, in the actuation formulation the pVl^ ] term is 

included since it is comprised of known quantities shown in Eq. (3.1 12). 

The actuation voltage used in the analysis to determine was applied to a single 

piezoceramic actuator attached to a 14"xl0"x0.040" clamped plate is shown in Figure 
5.12. The actuation signal used in the experiment was a 20- volt sine wave of 120 hertz 
obtained from a signal generator. The predicted and measured displacement of the plate 
are compared in Figure 5.13 and excellent agreement between the predicted displacement 
and test data is obtained. Since the displacement of test data is small, the signal-to-noise 
ratio of the accelerometer was less than ideal resulting in the excessive noise present on 
the displacement signal. 
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Wmax/h 



Fig. 5.1 Navier solution vs. finite element analysis non-dimensional displacement 
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Fig. 5.2 Single mode solution vs. finite element analysis non-dimensional displacement 





Fig. 5.4 Navier solution vs. finite element analysis sensor voltage 
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Fig. 5.6 Comparison of static sensor voltage from linear finite element analysis and 

experiment 
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Wmax/h 



Fig. 5.7 Comparison of non-dimensional displacement from nonlinear finite element and 

experiment 
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Fig. 5.8 Comparison of large deflection static sensor voltage from nonlinear finite 

element and experiment 
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Fig. 5.10 Comparison of dynamic sensor voltage output from finite element analysis and 

experiment 
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Fig. 5.12 Time history of electrical actuator signal used in the finite element analysis 



Finite Element Analysis 



Fig. 5.13 Comparison of actuator plate displacement from finite element and experiment 



CHAPTER VI 


CONCLUSIONS 

A finite element formulation exploiting the inherent electrical-structural coupled 
nature of piezoelectric materials is presented and compared to experimental test results. 
To verify the analytical model several tests were performed which included static uniform 
distributed loading, dynamic mechanical point loading, and dynamic piezoelectric 
actuation of a thin clamped plate with a bonded piezoceramic patch. The objective was to 
utilize linear piezoelectric theory in conjunction with the small deflection assumption in 
the finite element formulation, however the experimental tests indicated that the uniform 
distributed loading of the panel resulted in geometrical nonlinear transverse deflection. 
The predicted results obtained using small deflections were in excellent agreement with 
the test results for small pressure loading, however they were unacceptable for larger load 
values. The experimental pressure loading considered is not unreasonably high 
considering typical high speed flight vehicles experience pressures within the 3-5 psi 
range. Since the panel was selected to represent a typical aerodynamic structure, the 
finite element formulation was updated to include the von Karman large deflection theory 
which subsequently lead to coupled nonlinear electrical-structural stiffness. The 
predicted sensor voltages using the large deflection assumptions were in good agreement 
when compared to the experimental results. The sensor voltages obtained using the large 
deflection assumption were not as good as the predicted large deflection displacements, 
however the small deflection results are clearly unacceptable for the larger pressure 
loading. Excellent agreement between the predicted results and the dynamic test results 
were obtained utilizing the small displacement assumption since the experimental test 
used small dynamic loads. 

It was discovered that the electrical-structural coupling of the piezoceramic leads to 
coupled linear and nonlinear electrical-structural stiffnesses in the large deflection static 
sensor. The coupled stiffness appears in the Newton-Raphson tangent stiffness 
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formulation. Thus, even though linear piezoelectric theory was utilized the electrical- 
structural coupling was sustained through the large deflection geometric nonlinearity. 

Experimental measurement error was introduced since a charge attenuation circuit 
was required due to the large charge signal generated by the piezoceramic patch. The 
additional capacitance alters the low frequency response of the piezoceramic charge 
amplifier by decreasing the overall time constant thus, the measurement error may be 
reduced by modeling the effects of the charge attenuator during analysis. The charge 
attenuator may be eliminated from the experiment by utilizing a smaller piezoceramic 
sensor, thus compromising actuator performance, or utilizing a source follower to 
measure a proportional rate of strain in lieu of a proportional displacement measurement. 

A more accurate model of conductors and insulators would improve the electric field 
formulation in the large deflection static sensor finite element formulation. 
Improvements can be made to include more elements to accurately model the 
piezoceramic while utilizing a symmetrical finite element model. However, it was 
observed that it is important to properly model the electrical boundary conditions when 
using a symmetrical model. Since the static sensor experiment used a pressure load that 
approximated a step function, the static sensor results could be improved by simulating 
the actual transient pressure loading and computing the results through numerical 
integration. In addition a more accurate method of measuring the static plate deflection 
should be used in lieu of the dial indicator used. 

Future work may be applied in the area of incorporating piezoelectric thermal 
properties in the analysis or considering the application of control theory to the results. 
The piezoceramic physical properties show a strong temperature dependency which may 
be incorporated in the piezoceramic constitutive relations. Thus, a fully coupled 
electrical, structural, and thermal finite element model can be formulated. The finite 
element formulation was implemented in Matlab© which permits future opportunity to 
incorporate control and simulation analysis. 
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APPENDIX A 


TRANSFORMATION MATRICES 

The element used in the finite element model is a rectangular element consisting of 
twenty four structural degrees of freedom, accounting for bending and membrane 
displacements plus an additional electrical degree of freedom for each piezoelectric layer. 

The rectangular element consists of four comer nodes with displacements w, u, v and 
their derivatives w x , w y , w . The electrical d.o.f.'s are represented by the voltage V 

uniformly distributed over one side of the piezoelectric electrode, see Figure 3.1. 

Coordinate transformations are required to relate the local element coordinates to the 
global structural coordinates. The membrane nodal displacements, given by 

{ W „,} r =L U l “2 “3 “4 V 1 V 2 V 3 V ti ( A1 ) 

are obtained using the bilinear approximation functions in Eq. (3.6). The nodal 
membrane displacements can be obtained by substituting the element nodal coordinates 
into the appropriate approximation functions. Thus the membrane nodal displacements 
become 


' l (0,0,t) = b, 

(Al) 

'. 2 (a e ,0,t) = bi +b 2 a e 

(A2) 

• 2 (a e ,b e ,t) = b l +b 2 a e +b 3 b e +b t a e b e 

(A3) 

b e ,t) = b t +b 2 b e 

(A4) 

(0, 0,t) = b 5 

(A5) 

• 2 (a e ,0,t) = b 5 +b 6 a e 

(A6) 

• 3 (a,,b',t) = b 5 +b 6 a e +b 7 b e +b s a e b e 

(A7) 
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v t (0,b e ,t) = b 5 +b 3 b t (A8) 

where the element length and width are a e and b e respectively and the generalized 
coordinates are represented by the b coefficients. 

The membrane displacements in Eqs. (A1-A8) may be written in matrix notation as 
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Similarly the bending transformation is determined using sixteen d.o.f. Thus the sixteen 
bending nodal displacements are 

{w 6 f = {w, w 2 w 3 w 4 w, xl w, x2 w, x3 w, xi 

w,,i W„2 W >y3 W >»l W >,,2 (A10) 

The displacements are approximated using the cubic polynomial defined by Eq. (3.4). 
Thus, the derivative expressions are given by 

w, x ~a 2 + 2a 4 x+a 5 y + 3a 7 x 2 + 2 a s xy + a 9 y 2 +3 a u x 2 y 

+ 2 a n xy 2 +a l3 y 3 +3 a u x 2 y 2 +2 a 15 xy 3 +3 a 16 x 2 y 3 (All) 

w, y = a 3 + a 5 x + 2a 6 y + a s x 2 + 2a 9 xy + 3a XQ y 2 +a u x 3 

+ 2a n x 2 y +3a n xy 2 +2a u x 2 y +3a xs x 2 y 2 +3 a 16 x 2 y 2 (A12) 
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w, xy = a 5 + 2a z x + 2a 9 y + 3a u x 2 +4 a u xy 

+ 3 a u y 2 +6a H x 2 y + 6a l5 xy 2 +9a l6 x 2 y 2 +3 a l6 x 3 y 2 (A13) 

The nodal bending displacements are obtained by substituting the nodal coordinates into 

Eq. (3.4) and Eqs. (A4-A6). Thus the transverse displacements are given as 

w x (p s 0,t)^a x (A 14) 

w 2 (a e ,0,t) = a l +a 2 a e +a 4 a 2 + a 7 a] (A15) 

w,(a e ,b e ,t) = a t +a 2 a e +a 3 b e +a„a] +a 5 a e b e +a 6 b* 

+ a 9 a e b ] +a v) b] +a u a]b, +a n a)b] 

+ + a u a l b e + «i5«& 3 (A16) 

w i (0,b e ,t) = a 1 +a 3 b] +a 6 b, +a w b ’ (A17) 

Likewise, the slopes with respect to the x-axis are 

w, xl (0,0,t) = a 2 (A 18) 

w, x2 (a e ,0,t)-a 2 +2 a 4 a e +3 a 7 a 2 (A19) 

w, x3 {a e ,b e ,t) = a 2 +2 a 4 a e +a 5 b e +3 a 7 a] +2 a 8 a e b e 
+ a 9 b 2 +3 a u a 2 b e +2 a l2 a e b 12 +a u b 2 e 
+ 3 a u a 2 b 2 +2 a X5 a e b] + a X6 a]b] (A20) 

w, x4 ( 0,b e ,t) = a 2 +a 5 b e +a 9 b 2 +a n b 3 e (A21) 
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The slopes with respect to the y-axis are 

w, yl (0,0,t) = a 3 (All) 

w, y2 (a e >0,t) = a 3 +a 5 a e + a 8 a 2 + a u a 3 e (A23) 
w >y 3 (a e >b e> t) = a 3 +a 5 a e +2 a 6 b e + a 8 a 2 + la 9 a e b e 
+ 3 a l0 b 2 e +a n a] +la n a 2 b e + 3a u b 2 


+ 2 a u a]b e +3 a l5 a 2 b 2 +3 a 16 a 3 e b 2 (A24) 

(0,b e ,t) = a 3 +2 a 6 b e +3a l0 b 2 (A25) 

The rotations about the r-axis are 

w, xy i {0,0,t)-a 5 (A26) 

w, xy2 (a e >0,t) = a 5 +2 a 8 a e +3 a n a 2 (All) 

(a e ,b e ,t) = a 5 +2 a 8 a e +2 a 9 b e +3 a u a 2 

+ 4a u a e b e +3 a l3 b 2 +6a u e 2 b e 
+ 6a is a e b 12 + 9 a l6 a 2 b 2 (A28) 

w, xy4 (0,b e ,t)- a 5 +2 a 9 b e +3 a 13 b 2 (A19) 


Equation (A 14-22) may be expressed in matrix notion as {w b }~[T b ] 1 {a} where [T b ] _1 is 
defined as 
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APPENDIX B 

COORDINATE TRANSFORMATION 
B.l Transformed Reduced Stiffness Matrix 

The stress-strain relations for a thin orthotropic layer are of the form 


'12 J 


Qu Q\2 0 S \ 

Ql2 Q22 0 < S 2 > ( B1 ) 

. 0 0 1066 _ y 12 v 

where the subscripts l and 2 refer to material coordinate directions. The reduced stiffness 
matrix, \Q\ is a function of engineering constants as follows 


0ii = hl 


1-^12^21 


_ V \2^2 
i£12 — 


l“Vi 2 v 2 i 


0 




22 


1-^21 


Qa> ~ ^12 


(B2) 

(B3) 

(B4) 

(B5) 


The stresses in the material coordinate system can be transformed into the global x-y 
coordinate system by 



> __ 

cos 2 # 
sin 2 6 

sin 2 6 
cos 2 # 

-2sin0cos0 

2sin6cos0 

a 2 f 

Jxy. 


sinOcosO 

-sm&os# 

cos 2 #-sz>z 2 # 

J 12 . 


=K(e)\ 


\J*y 


(B6) 
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The strain transformation can be expressed as 


r > 

£ x 



> — 

y 


[r X y\ 



cos 2 9 
sin 2 9 

-2sin9cos9 


sin 2 9 
cos 2 9 
2sin9cos9 


sin9cos9 
-sin9cos9 
cos 2 9 -sin 2 9 


' £ \ ' 

< £ i > 

y 12 „ 


=[T c m 


1/3 J 


(B7) 


where 6 is the lamination fiber angle with positive rotation from the x-axis with respect 
to the principal material coordinate 7-axis as shown in Figure Bl. Combining Eqs. (Bl- 
B7) yields the laminate transformed reduced stiffness matrix 


where 


T 


xy J 


■fil. 

Q\2 

Qu~ 

£ x 

Ql2 

Q22 

Q26 

£ y 

_Q\(> 

Q26 

Q(>6 _ 

Yxy , 


q]=K(9)T\Q\tM 


(B8) 


(B9) 


B.2 Transformation of Piezoelectric f W" Constants 

Actuation strain shown in the stress-strain relations is proportional to the piezoelectric 
d constants and are presented in matrix form 




o 

0 

1^31 


0 

d, 


d 


^32 ^33 


15 

0 


(BIO) 


This research assumes thin piezoelectric layers polarized in the 3-direction and isotropic 
in the 7- and 2-directions. Thus the piezoelectric constant matrix is reduced to 


[d\ = 


d%\ dy2 


0 0" 
0 0 
0 


The actuation strain may be represented as 


(Bl 1) 
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(B12) 


{e\ = [d] r {E} = E 



Transformation of the principal material coordinates of the d constants to the global 
coordinates may be accomplished by utilizing the strain transformation 


[d 1 

X 



d 

y 

II 

d^2 

d 

l xy J 


0 


(B13) 


Thus for the isotropic assumption of the piezoceramic layers, 


fd 

X 


r j ' 

a 2l 

d 

y 

► = « 

d 2l 

d 

L xy J 


o 


(B14) 


103 




APPENDIX C 


ELECTROSTATICS 

C.l Introduction 

Electric charge establishes an electromagnetic field. The strength of the 
electromagnetic field is a function of the magnitude, location, velocity, and acceleration 
of the charges present. An electrostatic field results from stationary charges, hence is a 
function of position only. In the special case of electrostatics, note that the terminology is 
modified by dropping magnetic from the root electromagnetic. This infers that magnetic 
results from dynamic charge behavior. In essence, electric charges in motion (i.e. electric 
current) produce a magnetic field. 

C.2 Electrostatic Fields 

Coulomb's Law established experimentally, that a force exist between two charged 
bodies, which tend to repel or attract each other. The direction of the force depends on 
whether both bodies have similar or dissimilar charges. If two small spherical bodies, 
each of charge q j and q2 respectively, are present in an infinite homogeneous insulating 
medium separated by distance r, then Coulomb's force is expressed as 


4 7c€r 


(Cl) 


The direction of the force is along the line between the two charges, where e is called 
the permittivity or dielectric constant of the medium. In free space or a vacuum, £ is 
defined as £ 0 = 8 . 854 x 10 “^ Farad/meter, where the subscript o indicates free space. The 
permittivity of other materials are referenced to the permittivity of free space through 
their relative permittivity, that is £=e 0 £ r . In the MKS system, the unit of charge was 
previously defined as the coulomb (one ampere of current flowing for one second 
transports one coulomb of charge). Thus in order to leave Maxwell's field equations 


105 



independent of the constant 4tt, the units of Coulomb's force law were rationalized by a 
constant of 47T. 

In the above example imagine that the q 2 charge is fixed and consider qj as being 
available to move about, i.e. a test probe charge. As the test probe charge is located at 
various positions near the fixed charge it experiences a force. The magnitude and 
direction of the force depends on the location relative to the fixed charge. Hence 
surrounding the fixed charge there exists an electric force field. This electric force field 
is described to have a strength E, which has units of force per uni t charge. The 
magnitude of the force is given by Eq. (Cl) and the magnitude of the electric field 
strength is 


e =t ± t ( C 2 ) 

4 ner 

Note that the magnitude of the force is directly proportional to the magnitude of the test 
probe charge, however the electric field strength is defined as a force per unit charge thus 
independent of the magnitude of the test probe charge. Hence if the magnitude of the test 
charge approaches zero, the force diminishes; however the normalized force remains 
constant. Thus the electric field resulting from a charge exists, regardless of the presence 
of a test probe charge. In effect the test charge can be utilized to verify the existence of 
an electric field, and abstractly, if any charge q exists, then an electric field exists. The 
electric field around a point charge is a vector quantity and can be written as 

E = — rf (C3) 

4 ner 

where r is a unit radial outward vector from charge q. 
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C.3 Electric Displacement 

Equation (C3) indicates that the electric field is a function of the magnitude and 
position of the charge and the permittivity material. It is desirable to introduce the 
electric displacement or electric flux, that is independent of the material permittivity. The 
electric flux is best described by the results of Faraday's famous experiment with 
concentric spheres [16]. Faraday placed a metallic sphere with charge Q inside of another 
metallic sphere. Great care was taken to keep the spheres separated at all times. Then the 
outer sphere was grounded to earth for a very short time and then again with great 
diligence using insulating tools the inner sphere was removed. Once the outer sphere was 
reassembled, the charge on it was measured. The charge on the outer sphere was found to 
be equal and opposite in sign to the charge placed on the inner sphere. This experiment 
was repeated for several sphere sizes and with several dielectric materials and the results 
were always the same. Hence there is an electric flux or displacement through the 
dielectric. The displacement being a function of the magnitude of charge and 
independent of permittivity 


'¥ = Q 


(C4) 


Consider an isolated point charge far from any other particle. The electric displacement 
density or electric displacement per unit area on any point of a sphere surrounding the 
point charge is 


D = 


¥ 

Am 2 


_JL 

Am- 2 


(C5) 


The displacement density is a vector quantity with its direction outward normal to the 
sphere 


D = 


¥ - 

T r 

Am 


(C6) 
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Upon examination of Eq. (A.3) and Eq. (A.6), D and E are related and can be written as 
the electromagnetic constitutive relation 


D = eE (C7) 

where the permittivity is a tensor for the most general anisotropic materials 



C,i 

^12 

fc 13 

€ = 

^12 

^22 

£ 23 


_ e i3 

^23 

^33 _ 


(C8) 


C.4 Potential Function 

The electric field is a conservative force field. Hence there is no dissipation 
mechanism and all energy must be stored in either potential or kinetic form. The work 
done, on the system above, in moving the test probe charge around the fixed charge, i.e. 
against the force F, can be calculated as 


work = - £ Fdr (C9) 

As with any conservative force field an arbitrary reference may be chosen. As shown 
above, infinity is commonly selected as the reference. Thus the work done on the test 
probe charge is 


work = ~ f d r 
Aner 


<h<h 

AneR 


(CIO) 


The normalized work done on the test probe charge is defined as the electric potential at a 
point due to the presence of charge q j, hence 


V = 


work 

unit charge 


AneR 


(Cll) 
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Note that potential is a scalar quantity. The potential has a unit of joule per coulomb 
which is the MKS units for a volt. For electrostatics, the terms potential, potential 
difference and voltage are synonymous. Another useful relationship is the differential 
work or differential voltage required in moving a positive unit charge an infinitesimal 
distance. For example, the differential work is 


dW = dV ~-E‘ds 


(C12) 


where dV can be expressed as VV • ds , thus E = -VV and the electric field strength at any 
point is simply the negative of the potential gradient at the point. If Eq. (Cl 2) is 
integrated the potential difference will become much clearer. For example, consider the 
potential difference between two points a and b as 

= (C13) 

Thus the potential difference between point a and b is 

V a -V b = \ *E-<* (Cl 4) 

C.5 Capacitance 

The capacitance between two conductors is defined as the ratio of charge on the 
conductors, to the voltage or potential difference between each of the conductors 



(Cl 5) 


Note that the capacitance is defined irrespective of the size and shape or distance between 
the conductors. Gauss's law can be used to verify that the total charge on one of the 
conductors is indeed proportional to the potential difference, and the proportionality is 
capacitance. 
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C.6 Electrostatic Energy 

When a capacitor is charged, a voltage or potential difference exists between the 
conductors. As previously discussed, the voltage will establish a conservative electric 
force field between the two conductors. Thus energy storage is present in a charged 
capacitor, since no dissipation mechanism exists within a conservative force field. The 
amount of energy stored, is simply the work required to charge the capacitor. Recall that 
potential was defined as work done per unit charge in Eq.(Cl 1), thus the differential work 
is 

dW = Vdq (Cl 6) 

and the total work is computed as 



The stored energy of a charge capacitor can also be written as 


(Cl 7) 


U = -V 2 C (Cl 8) 

2 

Another approach to visualize the energy stored or work done in charging a capacitor is to 
consider the necessary energy required in establishing a charge distribution in space. If a 
free space is considered, and N discrete charges are brought into a given volume, then the 
work done or energy expended in locating the charges will be 


Q Q 


1 


w M R 0 


i*j 


(Cl 9) 


If a continuous charge distribution is considered instead of discrete charges, Eq. (Cl 9) 
can be written in integral form as 


W = 


1 

8 7T€ 


n 


pM)Pc^i) 

R 


d¥ x d,V . j 


(C20) 
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where qf and qj were replaced with volume charge densities p c (r x )d¥- x and p c (r 2 )dV- 2 
respectively, and R = \r x - r 2 1 . Recall that potential is defined as work per unit charge in 
Eq. (Cl 1), so work done on the charge distribution can be written as the potential by 


V(r) = 


W_ 

/>>,) 


= — (C21) 

4 Tie J R 2 

Substituting Eq. (C21) into Eq. (C20) expresses the work necessary to create a continuous 
charge distribution 


\p{ r y(r)dV (C22) 

If the vector identity V • (FD) = W D + D VV is applied to Eq. (C22) it may be re- 
written as 


W = - \pVdV- 

2 fc 


= ~W ■ OdV- 
2 


= - J[V(FD)-D-VFj/^ 

2 r z 


-[VD-dn + -\D-'Edi L (C23) 

2 s ^ v- 


For a piezoelectric generator, the externally applied voltage is zero and a mechanical 
deformation or strain results in a Maxwell's self induced electric field. The electrical 
energy density thus becomes 

= (C24) 

2 J 

^ v- 
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This work formulation assumes that there is no free or space charge present within the 
piezoceramic. 
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APPENDIX D 


PIEZOCERAMIC ADHESIVE 

The adhesive used to bond the piezoceramic patch to the panel was Micro- 
Measurements GA-2 strain gage adhesive system [23]. The GA-2 system is an epoxy 
adhesive that is combined with the hardener 10-A. The adhesive is specified to have 10% 
to 15% elongation capabilities after 40 hours of curing time at 21° C. It is recommended 
that a clamping force of 5-20 psi be applied to the patch during the curing process to 
assure a complete bond. Once the epoxy is mixed with the hardener the pot life is 
approximately 15 minutes at 21° C. The adhesive is a dielectric which prohibits 
electrical conduction between the piezoceramic and the plate. 
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APPENDIX E 


CHARGE AMPLIFIER DATA 

The charge amplifier used, model number 422M77, was provided by PCB 
Piezotronics, Inc. [27]. The gain is fixed at 0.2537 mV/pC with a bias voltage of 11.09 
V, a feedback capacitance of 4000 pF and a feedback resistance of 3x1 0^ ohms. 
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